Hall Transport in Granular Metals and Effects of Coulomb Interactions 
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We present a theory of Hall eflFect in granular systems at large tunneling conductance 3> 1- 
Hall transport is essentially determined by the intragrain electron dynamics, which, as we find using 
the Kubo formula and diagrammatic technique, can be described by nonzero diffusion modes inside 
the grains. We show that in the absence of Coulomb interaction the Hall resistivity p^j^ depends 
neither on the tunneling conductance nor on the intragrain disorder and is given by the classical 
formula pxy = H/{n*ec), where n* differs from the carrier density n inside the grains by a numerical 
coefficient determined by the shape of the grains and type of granular lattice. Further, we study 
the effects of Coulomb interactions by calculating first-order in l/gr corrections and find that (i) in 
a wide range of temperatures T > F exceeding the tunneling escape rate F, the Hall resistivity pxy 
and conductivity a^y acquire logarithmic in T corrections, which are of local origin and absent in 
homogeneously disordered metals; (ii) large-scale "Altshuler-Aronov" correction to Oxy, relevant at 
T ^ F, vanishes in agreement with the theory of homogeneously disordered metals. 

PACS numbers; 73.63.-b, 73.23.Hk, 61.46.Df 



I. INTRODUCTION 

Hall transport in different systems has always been 
a subject of extensive research. Already the classical 
Drude-Boltzmann theory provides us with an interesting 
result. It is well-known that the Hall resistivity (HR) 
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(1.1) 



of a disordered metal does not depend on the mean free 
path and is determined solely by the carrier concentra- 
tion n allowing one to extract it experimentally. At low 
enough temperatures quantum effects of Coulomb inter- 
action and weak localization (see, e.g., Refs. 1,2) influ- 
ence the Hall transport, giving corrections to Eq. (|1.1[) . 

Dense-packed arrays of metallic or semiconducting 
nanoparticles imbedded into an insulating matrix, usu- 
ally called granular systems or nanocrystals, have re- 
cently received much attention from both experimen- 
tal and theoretical sides (see a Review ^ and references 
therein). The longitudinal transport in such systems is 
theoretically well understood now, both in the metallic 
and insulating regimes. At the same time. Hall trans- 
port in such granular materials has not been addressed 
theoretically before, neither in the insulating nor in the 
metallic regimes. The absence of a theoretical descrip- 
tion is apparently one of the reasons, why measurements 
of the Hall resistivity have not become a standard tool 
for characterization of granular metals, although they do 
not seem to be very difficult. 

Trying to apply the conventional theory of disordered 
metals to granular systems, the following questions can 
be asked: To what extent is the formula (|l.ip applica- 
ble to granular metals? How is the carrier concentration 
extracted from Eq. (|1.1[) related to the actual carrier con- 
centration inside the grains? What impact do quantum 
effects have on Hall transport of a granular system? 
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FIG. 1: Granular system and classical picture of Hall con- 
ductivity. The external Ohmic voltage Vy is applied to the 
contacts in the y direction. The resulting Ohmic current 
ly = GrVy running through the grain in the y direction causes 
the Hall voltage drop Vh = Rnly between its opposite banks 
in the x direction. Since for calculating Hall conductivity o^y 
the total voltage drop per lattice period in the x direction is 
assumed 0, the Hall voltage Vh is also applied (with an oppo- 
site sign) to the contacts in the x direction, causing the Hall 
current I^ = GtVh = GlRnVy [Eq. (fT31l ]. 



In this paper we present a theory of Hall effect in a 
granular system in the metallic regime and answer these 
questions. 

In the metallic regime ( granular metaF ) , when the in- 
tergrain tunneling conductance Gt = (2e^/?i)(7T is large, 
(7t ^ 1 (further we set h = 1), the granular system as 
a whole is, roughly speaking, a good conductor and its 
properties are quite similar to those of ordinary homo- 
geneously disordered metals (HDMs). At the same time 
the granularity of the system brings a new physical as- 
pect: confinement of electrons inside the grains. In a 
system with "well-pronounced" granularity electron tra- 
verses each grain many times before it escapes from it to 
some neighboring grain due to tunneling. This is ensured 
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by the condition that the tunnehng escape rate F is much 
smaher than the Thouless energy _Btii of the grain: 

r«£;Th, (1.2) 

or, equivalently, the tunnehng conductance Gt is 
much smaller than the longitudinal conductance Go — 
{2e'^/h)go of the grain: 

ffT<5o, (1-3) 

since T — gxS and £^Th oc god {d is the mean level spacing 
of the grain) . 

The conditions p.2|) and ()1.3|) , leading to new physic^ 
absent in HDMs, simphfy calculations at the same time. 
Consider, for example, the classical (in the absence of 
quantum effects, such as Coulomb interaction and weak 
locahzation) longitudinal conductivity (LC) a^J of a reg- 
ular quadratic or cubic granular lattice (Fig. [T|) with all 
contacts having equal conductances Gt- In the limit 
gr ^ 9o the main contribution to the longitudinal re- 
sistivity (LR) p^xx = l/cra comes from the tunnel bar- 
riers between the grains rather than from scattering on 
impurities inside the grains and LC equals 

aW^Gra^-^ (1.4) 

where a is the size of the grains and d = 2, 3 is the dimen- 
sionality of the array. The longitudinal conductance Go 
of the grain itself, which, in principle, should be obtained 
from a solution of a classical electrodynamics problem for 
the distribution of the electric potential inside the grain 
and is, therefore, determined by the properties of the in- 
tragrain electron dynamics, does not enter Eq. (II. 4p . 

Thus, when studying longitudinal transport one may 
neglect the details of electron dynamics inside the grains, 
which is a significant simphfication. Technically, this 
is equivalent to considering only the zero (coordinate- 
independent inside the grains) spatial modes of the dif- 
fusons or phases in the phase functional^i^. Owing to the 
conditions (II. 2p and (jl.Sp . the zero-mode approximation 
suffices for studying the longitudinal transport. 

For Hall transport, however, the situation appears to 
be more comphcated. The Hall current originates from 
the transversal drift in the crossed magnetic and elec- 
tric fields inside the grains. The classically prohibited 
regions of tunnel contacts are neglegibly small for dense- 
packed arrays and cannot contribute to the Hall trans- 
port. From simple classical considerations (see Fig. [T]) 

one obtains, that Hall conductivity (HC) a^J in the lead- 
ing in gx/go 1 order is 

ai") = G|i?ffa2-^ (1.5) 

where Rh is the Hall resistance of the grain. Just like 
Go, the Hall resistance Rh should be obtained from a 
solution of a classical electrodynamics problem for the 
distribution of the electric potential inside the grain. We 
come to the situation when one is forced to take the in- 
tragrain electron dynamics into account, no matter how 



well the conditions (|1.2p and (|1.3p are satisfied. In other 
words, the zero-mode approximation is not sufficient for 
the description of the Hall transport of a granular system. 

However, a purely classical approach to the problem, 
giving a quick answer (|1.5p , does not allow one to include 
quantum effects (such as Coulomb interaction and weak 
locahzation) into considerations, which come into play at 
sufficiently low temperatures and can significantly affect 
transport properties. 

In this work we develop a method of calculating con- 
ductivity of a granular system in the metallic regime, 
which allows one to take the intragrain electron dynam- 
ics into account. Using the Kubo formula and diagram- 
matic technique, we show that this can be done by consid- 
ering nonzero {coordinate- dependent) modes of standard 
two-particle propagators ( "diffusons") inside the grains. 
This procedures accounts for the finiteness of the ratio 
9t / go and reproduces the solution of the classical elec- 
trodynamics problem for the conductivity of a granular 
medium. The generality of our approach allows one, in 
principle, to study both LC and HC of the granular sys- 
tem for arbitrary ratio gr / go and for arbitrary type of 
the intragrain electron dynamics, either ballistic or dif- 
fusive. Nonzero modes of the diffusons are eventually 
related to the longitudinal Gq^ and Hall Rh resistances 
of the grain. 

We apply our method to the problem of Hall trans- 
port, for which considering intragrain dynamics is in- 
evitable. Neglecting quantum effects, we do recover the 
classical formula (II. 5p for the Hall conductivity and ob- 
tain quite a universal result for the Hall resistivity. Dia- 
grammatic approach aUows us to include quantum effects 
of Coulomb interaction and weak localization straight- 
forwardly into the developed scheme. We study the in- 
fiuence of Coulomb interactions on HC and HR by cal- 
culating first-order corrections and find that the major 
temperature dependence of both HC and HR of a granu- 
lar metal comes from the contributions which are absent 
in HDMs. We also announce our results for weak local- 
ization corrections, detailed calculations of which wiU be 
presented elsewherei^. Part of the results of our work 
(for temperatures T > T) was presented in a brief form 
in Ref . % 

The paper is organized as foUows. In Sec. |TT] we 
present our results for Hall conductivity and resistivity 
and Coulomb interaction corrections to them. In Sec. IIIII 
the model for the granular system is formulated and dis- 
cussed. In Sec.|IV]the main features of the diagrammatic 
technique are explained, and important building blocks, 
namely, the intragrain diffuson in the presence of mag- 
netic field and the screened Coulomb interaction, are ob- 
tained. In Sec.|V]we calculate Hall conductivity neglect- 
ing quantum effects of Coulomb interaction and obtain 
the correspondence with the classical result. Quantum 
effects of Coulomb interaction are studied in Sec. IVII 
Concluding remarks are presented in Sec. IVIII In Ap- 
pendix [X] the boundary condition for the intragrain dif- 
fuson in the presence of magnetic field is derived. 
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II. RESULTS 

In this section we list the main results of this work. 
We perform calculations for magnetic fields H such that 
'^hTq <C 1, where ujh = eH/(mc) is the cyclotron fre- 
quency and To is the electron scattering time inside the 
grain. Since the (effective) electron mean free path 
I = vpTo < a does not exceed the grain size a, and 
typically a « 1 — lOOnm, the condition lohtq ^ 1 is 
well fulfilled even for experimentally high fields H. We 
also assume that the granularity of the system is "well- 
pronounced", i.e., the conditions ()1.2p and (jl.3p are sat- 
isfied. Other assumptions and approximations are for- 
mulated in Sec. IIIII 

Classical Hall conductivity and resistivity. First, we 
neglect quantum effects of Coulomb interaction and ob- 
tain Eq. (jl.5[) for HC a^J in the lowest nonvanishing or- 
der in gx/go ^ 1. This result obtained by diagrammatic 
methods is of completely classical origin provided the 
tunneling contacts are viewed as surface resistors with 
conductance Gt- The HR of the system 
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(2.1) 



following from Eqs. (jl.4p and p.5p . thus, does not depend 
on the tunneling conductance Gt and is expressed solely 
through the Hall resistance Rh of a single grain. Further, 
the Hall resistance Rh of the grain does not depend on the 
intragrain disorder, but only on the geometry of the grain 
and carrier density n of the grain material. For grains of a 
simple geometry (e.g., having reflectional symmetry in all 
three dimensions) Rh = P^yO,/S , where p|y = H/{nec) 
is the specific Hall resistivity of the grain material and S 
is the area of the largest cross section of the grain. 

Therefore, akin to the universal result (jl.ip for ordi- 
nary disordered metals, for the classical Hall resistivity 
of a granular metal we obtain 



p(0) 
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(2.2) 



in the case of a three-dimensional sample (3D, d ~ 
many grain monolayers). Here, 

S 

n* = An, v4 = — < 1, 

is the effective carrier density of the system, which dif- 
fers from the actual carrier density n inside the grains 
only by a numerical factor A determined by the shape of 
the grains [A = 7r/4 for spherical and A = 1 for cubic 
grains). For a two-dimensional sample (2D, d = 2, one 
or a few grain monolayers) the expression (|2.2p must di- 
vided by the thickness dz of the sample or, equivalently, 
n* — dzAn in this case'^^. 

The result (|2.2p for the Hall resistivity pS] is quite uni- 
versal. It is valid even if (i) the tunneling conductances 
Gt fluctuate from contact to contact and (ii) the mean 



free path I fluctuates from grain to grain: HR is sim- 
ply independent of the distributions of Gt and I; There- 
fore, Eq. (|2.2p is applicable to real granular arrays in 
which such irregularities are always present (provided 
such system is still in the metallic regime). We also 
note that although Eq. (|2.ip was obtained for a regular 
quadratic/cubic granular lattice, the result (|2.2p with a 
different factor A < 1 remains valid for other regular lat- 
tices (e.g., more common for real experimental samples 
triangular lattice). We also expect Eq. (|2.2p to hold for 
arrays with moderate structural disorder, i.e., in which 
the positions of the grains deviate from regular and their 
sizes and shapes are not identical. 

Coulomb interaction corrections. Next, we calculate 
the first-order in l/gr corrections to HC a^J [Eq. p.Sp ] 
due to Coulomb interaction. We find significant contri- 
butions for temperatures T < gTEc not exceeding the 
inverse RC time gTEc of the system "grain-|-contact" 
[Ec — I (na) is the charging energy and k is the dielec- 
tric constant of the array], whereas for T > gTEc the 
relative corrections are of the order of l/gT or smaller. 
Three types of corrections to HC a^xj [Eq. (|1.5p ] can be 
identified: 



.(0) 



■6a^y^ + Sa^^^ + Sa^f. (2.3) 



The first correction Sa'^^ can be attributed to the renor- 
malization of the individual tunneling conductances Gt 
[tunneling anomaly (TA)^'^'^] in the granular medium 
and has the form 



^ xy 
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TigTd 



;ln 



grEc 



max(T, F) 



for T < gTEc 



(2.4) 

This correction renormalizes the tunneling conductances 
Gt in Eq. (|1.5p . but does not affect the Hall resistance 
Rh of the grain. 

The second correction Sa^^ corresponds to the process 
of virtual diffusion (VD) of electrons through the grain 
and equals 
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In 



mm{gTEc, ^Th) 
max(r, F) 



(2.5) 



for T < min((7T-E'c, i^Th), where ~ 1 is a numerical lat- 
tice structure factor (|6.23p . Contrary to Scr"^^, the cor- 
rection Sa^^ is suppressed at temperatures greater than 
the Thouless energy of the grain i?Th, which emphasizes 
its diffusion character. 

For T > F both corrections ba^y and ba^^ are 
InT-dependent. This dependence saturates at temper- 
atures r ~ F, so that (Jtrjj^ and Scr^^ remain logarith- 
mically large constants at T < F. These two correc- 
tions are specific for granular systems and, in essence, 
due to the strong discrepancy of timescales of the intra- 
(^Th) ^^'^ intergrain (F~^) electron dynamics described 
by Eq. (|1.2p . They arise from spatial scales of the or- 
der of the grain size a and are absent in HDMs. The 
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logarithmic behavior of the corrections (|2.4|) and (|2.5p is 
due to the form of the screened Coufomb interaction in 
granular systems^. They have the same logarithmic form 
in 2D and 3D, but the coefRcients are not universal and 
lattice-dependent: 1/d and Cd [Eq. (|6.23[) ] are the results 
for the cubic (3D) or quadratic (2D) lattice, which we 
assumed in our calculations. 

The third correction 5(j-^.^ in Eq. (|2.3p is analogous 
to the one present in homogeneously disordered metalsi 
[ "Altshuler-Aronov" ( AA) corrections] . It can be signifi- 
cant at low temperatures T <C F only, when the thermal 
length Lt = \/ D^/T ^ a for the intergrain motion ex- 
ceeds the grain size a {Dq = Ta^ is the effective diffusion 
coefficient for the intergrain electron motion at scales 
greater than a). However, we find that this correction 
vanishes identically both in 2D and 3D: 



da. 



AA 
xy 



0. 



(2.6) 



It is always instructive to compare the results for a gran- 
ular metal with those for a HDM. The quantities arising 
from large spatial scales (exceeding the grain size a for 
a granular metal and the mean free path I for a HDM) 
are expected to behave universally, because at such scales 
the microscopic structure of the system becomes irrele- 
vant. Indeed, the result (|2.6p for Scr^,^ agrees with the 
one for HDMs first obtained in Ref. 8. Being an exact 
cancellation, however, Eq. (|2.6p is valid not only for low 
T r, but for arbitrary relevant temperatures. 

The quantity directly measured in experiments is the 
Hall resistivity 



3D case, being smaller than the logarithmic contributions 
(|2.4p . (|2.5p . (|2.9p for all relevant temperatures down to 

very low ones^^: 5a^^/ai°] < I /gr- 
in 2D case^^, i.e. for granular films of thickness dz 
consisting of one or a few grain monolayers {dz/ a is the 
number of monolayers), the correction 5a^^ is diverging 
at large spatial scales. This divergence is relevant for low 
temperatures T ^ r(a/c?2)^ (when ^ dz), for which 
Sa^^ acquires a logarithmic dependencoi^: 



^ XX 
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r / a 

T Uz 



The total Coulomb interaction correction 

„TA , xVD , x„AA 



(2.10) 



(2.11) 



to Hall resistivity (|2.7p is given by the sum of the con- 
tributions arising from the corresponding corrections to 
Hall (|2.3p and longitudinal (|2.8p conductivities as 

dpxy _ daxy ^ TA,VD,AA. 



pxy ^ 



^ XX 



Since the TA effects leads to the renormalization of the 
tunneling conductance Gt only, it cannot affect the HR 
p^xy [Eq. (|2.2p ]. which does not contain Gt- Indeed, it 
follows from Eqs. and ([^ that 



TA 



XX 



^xy (Q) , c- 

P^V^Zr= Pxy +S Pxy, 



(2.7) 



where pl^/j is the "bare" HR [Eq. (E^])], Sp^y is the to- 
tal Coulomb interaction correction to HR and HC axy 
is given by Eq. (|2.3p . The interaction corrections to LC 
were studied in Refs. |l[l3 and the following result was 
obtained: 



r TA r 
OCT,r^ -I- 0(7, 



AA 



(2.8) 



[(5ctJ/ and 5a^^ correspond to Sai [Eq. (2b)] and Sa2 
[Eq.(2c)] in Ref. [l^, respectively]. The correction (5crJ/ 
is due to the tunneling anomaly in granular metal. It 
renormalizes the tunneling conductance Gt in Eq. (|1.4p 
and equals 



XX 



2TTgTd 



:ln 



9tEc 



max(T', r) 



for T < gTEc- 



(2.9) 

This correction is of local origin and governs the temper- 
ature dependence of LC axx{T) in a wide temperature 
range. The Hall counterpart of 5crJ/ is Sa'^,^ [Eq. 1^^ ]. 

The correction Sa^^ is analogous to that in a HDM, 
first obtained by Altshuler and Aronov (AA) in Ref. [ill 
and its HaU counterpart is 6a^^ [Eq. ^^]. The AA 
correction Sa^^ does not diverge at large spatial scales in 



and, therefore, the correction to Hall resistivity from the 
"tunneling anomaly" effect vanishes: 

SPxy = 0. (2.12) 

Further, since the "virtual diffusion" correction is absent 
for longitudinal conductivity [Eq. (|2.8p }^ we obtain 



Pxy 



,(0) 



(2.13) 



where Sa^^{T) is given by Eq. (|2.5p . Finally, since the 
"Altshuler-Aronov" -type correction Jcr^jJ* to Hall con- 
ductivity vanishes [Eq. (|2.6p ] , we get 



(2.14) 



Pxy 



= -2- 



XX 



Therefore, according to Eqs. (PTT|) . ((TT^ . ([^T^ . and 
p.l4p . the total Coulomb interaction correction Sp^y to 
Hall resistivity is 



SpxyjT) Sa^y^jT) JatHT) 



Pxy 



(Jxy 



(y XX 



(2.15) 



where Sa^y" is given by Eq. (23]) and Sa^^ [Eq. ((TTU)) ] 
can be significant in 2D at T ^ T{a/dz)'^ only. 
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FIG. 2: Temperature dependence of the Coulomb interaction 
corrections 5ply{T) [Eq. fTH^ ] and 5pi^{T) [Eq. d^T^] t o 
the Hall resistivity [second and third terms in Eq. (|2.17|l . 
respectively]. The most significant contribution in a wide 
range of temperatures [T < min((/T£^c, -Exh)] both for two- 
and three-dimensional samples comes from the correction 
Sp^^iT) (thick solid line), which is due to the process of 
"virtual diffusion" of electrons through a single grain. The 
contribution SpY:^ (T) is of local origin and absent in homo- 
geneously disordered metals. It depends logarithmically on 
temperature T in the range F < T < min((;T-Ec, -Etii)) sat- 
urating at T ~ F and remaining constant for T < F. The 
correction &pxy'{T) is analogous to the "Altshuler-Aronov" 
correction in ordinary disordered metals, it can be significant 
[Eq. (|2.10|l ] for sufficiently thin granular films and low enough 
temperatures [T ^ Y'{a/dz)'^\ only, the latter 2D case shown 
in figure (dashed line). 



Equations ([^ -(1 ^ . (PTTjl - fTTS)) constitute our 
main result for Coulomb interaction corrections to Hall 
conductivity and resistivity. Another effect occurring 
at similar temperatures is weak localization (WL). The 
WL corrections to LC of a granular metal were stud- 
ied in Refs. IT^fTslfl^ . Significant (logarithmic) contribu- 
tions may arise in 2D samples only, from spatial scales 
greater than the grain size a, when the inverse dephas- 
ing time < F (if l/r^ cx T/gy ^^d^ , this corresponds 
to T < gr^)- However, we findi^ that the first-order 
in 1/(?T WL correction to Hall resistivity vanishes iden- 
tically both in 2D and 3D in correspondence with the 
result for HDMs^ii^JI: 



(2.16) 



Therefore weak localization does not affect Hall resistiv- 
ity at least in the first order in l/gxi in which we obtain 
significant contributions from the Coulomb interaction. 

Summarizing our findings for the classical result (j2.2p . 
Coulomb interaction [Eqs. (HHl-dMl) (|2lT|) - ((2d5)) ] and 
weak localization (|2.16p corrections, we predict the fol- 
lowing behavior of the Hall resistivity of a granular metal: 



PxyiT) 
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min (ffT-^c, Etii) 
max(T, F) 



^pZ' = 



^ XX 



(i) At high enough temperatures T > min(gT-E'c, ^'xh) 
the Hall resistivity pxy{T) = p^y is given by the Drude- 
like expression (|2.2p [first term in Eq. (|2.17p ] and is in- 
dependent of both the intragrain and tunnel contact dis- 
order. Measuring pxy at such T and using Eq. (|2.2p one 
can extract an important characteristic of the granular 
system: its effective carrier density n* = An, which is re- 
lated to the actual carrier density n of the grain material 
through a geometry-dependent factor ^ < 1. 

(ii) In a wide temperature range E < T < 
mm{gTEc, -ETh) both for 2D and 3D samples, local effects 
of Coulomb interaction lead to the logarithmic in T cor- 
rection to the Hall resistivity [5pY^{T), second term in 
Eq. (PT7|) . see Eqs. This In T-dependence 
saturates at T ~ F and 6pY^{T) remains constant for 
2^ ^ r. We emphasize that this correction is absent in 
homogeneously disordered metals, but it appears to be 
the major quantum correction to the Hall resistivity of a 
granular metal that governs the T-dependence of pxy {T) 
in a wide temperature range both for 2D and 3D samples. 

(iii) An additional T-dependence of pxy{T) may arise 
due to the "Altshuler-Aronov" correction 5(7^^{T) to 
the longitudinal conductivity [5p{^^{T), third term in 



(2.17) 



Eq. ((2T7ll . see Eqs. ([2:6)) . ((2d0| . (12141 ] at much lower 
temperatures T ^ F, the most significant logarithmic 
contribution expected for sufficiently thin granular films 
[one or a few grain monolayers, Eq. (|2.10p ]. 

The temperature behavior of the contributions 
Sp^y^iT) and Sp^^iT) is shown in Fig. [3 

We expect our result Eq. (|2.17p to hold for realistic 
arrays with moderate structural disorder and, most im- 
portantly, with randomly distributed tunneling conduc- 
tances, which is inevitable in real systems: (i) n* simply 
does not depend on the distribution of Gt] (ii) the loga- 
rithmic form of the major quantum correction [SpY^{T), 
second term in Eq. (|2.17p ] persists in this case, although 
the structure factor Cd ^ 1 does depend on the distri- 
bution of conductances and gx should be substituted by 
some averaged quantity. 

Comparison of our findings with experimental data 
may serve as a good check of the theory developed here. 
The experimental situation related to our theory is men- 
tioned in the Conclusion. 



III. MODEL AND HAMILTONIAN 

We consider a quadratic {d — 2, 2D) or cubic {d = 3, 
3D) lattice of metallic grains coupled to each other by 
tunnel contacts (Fig. [1]). 

Aiming to concentrate on the Hall effect, we assume 
the simplest case of translationally invariant lattice, i.e., 
equal tunneling conductances Gt of all contacts, transla- 
tionally invariant capacitance matrix, and identical prop- 
erties of all grains (the same form and size, mean free 
path, electron density, density of states, etc.). After the 
main properties of the Hall effect in such system have 
been established, we argue that they also hold for realistic 
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arrays. In real systems in the metallic regime, the major 
type of irregularities that (could) affect electron trans- 
port even for structurally quite regular arrays seems to 
be the randomness of tunneling conductances Gt, while 
other assumptions can be well met or are inessential. 

To provide more explicit analysis we further simplify 
the calculations assuming the intragrain electron dynam- 
ics diffusive, i.e., that the bulk mean free path / in the 
grain is much smaller than the size a of the grain, I <C a. 
In this case details of electron scattering off the grain 
boundary are irrelevant. However, our approach is also 
entirely applicable to the case of ballistic {I > a) intra- 
grain disorder, when surface scattering becomes impor- 
tant. The main results, listed in Sec. HH are valid for 
both diffusive and ballistic grains. 

In the metallic regime {gr S> 1) quantum effects of 
Coulomb interaction can be considered perturbatively 
with an expansion parameter I/qt as long as the rela- 
tive corrections remain small. 

We write the Hamiltonian describing the system as 



H — Ho + Ht + He- 
In Eq. ((3?T1) the first term 

^0 = E / *iV'^(ri) (pi - ^A(ri 



(3.1) 



(3.2) 

is the Hamiltonian of isolated grains, ^(p) = p^/(2m) — 
ep, ■A-(ri) is the vector potential describing the uni- 
form magnetic field H — Hgz directed along the z axis, 
U{Yi) is the random disorder potential of the grains, 
i = (ii, . . . , id) g Z'' is an integer vector numerating 
the grains. The integration with respect to ri is done 
over the volume of the grain i. Since we do not deal 
with spin-related phenomena in this paper, we omit the 
spin indices of the operators il^{v\). Accounting for spin 
degeneracy in the course of calculations is simple: each 
electron loop comes with the factor 2. We consider white- 
noise disorder and perform averaging using the Gaussian 
distribution with the variance 



{U{v,)U{v[))u 



1 



2'kvtq 



Sin - rO, 



(3.3) 



where v is the density of states in the grain at the Fermi 
level per one spin projection and tq is the scattering time. 
The tunneling Hamiltonian Ht in Eq. p.ip is given by 



H 



= ^(Xi,j-KXj,i), 



(3.4) 



(iJ) 



where the operator Xi j describes tunneling from the 
grain j to the grain i, the summation is taken over the 
neighboring grains connected by a tunnel contact, such 
that each contact is counted only once. For studying Hall 
effect the geometry of the grains and contacts is essen- 
tial, therefore we write the tunneling operators Xij in 
the coordinate representation: 



y cfsidsj i(si,Sj)i/''^(si)V'(sj), 



(3.5) 



where the integration is carried out over two surfaces 
of the contact: one of them (si) belonging to the i-th 
grain, whereas the other (sj) to the j-th grain. Such 
form implies that tunneling occurs from a close vicinity 
of the contact of atomic size, but not from the bulk of the 
grain. This is a natural assumption, since we consider a 
good metallic limit for the grains, i.e., the size a of the 
grains is much greater than the Fermi length, ppa ^ 1 
{pf is the Fermi momentum). Fast oscillations of the 
wave functions in the grains result in a fast decay of the 
overlap of the wave functions in different grains. Since 
h} — Ht, we have X^^ — Xyi and t*{si, sj) = i(sj, Si). 

Without further assumptions about the tunneling am- 
plitudes i(si,Sj) in Eq. (|3.5p . electrons can tunnel from a 
given point sj to an arbitrary point Si on the other side of 
the contact. However, it is physically clear that (i) elec- 
trons effectively tunnel from the point sj to the points Si 
in the vicinity of sj of atomic size only, therefore t{si, sj) 
should decay rapidly on atomic scale as a function of 
Si — sj; (ii) the amplitude i(si, sj) can also fluctuate as a 
function of Si for fixed Si — sj due to irregularities of the 
contact on atomic scale. 

To effectively model this behavior of the tunneling am- 
plitudes we consider i(si,Sj) as Gaussian random vari- 
ables and average over them with the variance 



(i(si,Sj)i(sj,Si))t = tlS{si - Sj), 



(3.6) 



where S{si — sj) is an atomic scale (5-function on the con- 
tact surface and has a meaning of tunneling probability 
per unit area of the contact. As we will see, the assump- 
tion ppa ^ 1 will enable us to neglect the contributions 
containing the regular parts (i(si,Sj))t of the timneling 
amplitudes. 

The third term in Eq. (|3.ip stands for the Coulomb 
interaction between the electrons. In principle, one has 
to start with its the bare form 



Hr 



-E 

2 ^ 

iJ 



dridvj i/''^(ri)V'^(i"j)i 



-?/'(rj)V'(ri). 



(3.7) 

Proceeding with the calculations we will have to take the 
screening of Coulomb potential by electron motion into 
account. One should distinguish between the intragrain 
and intergrain electron motion. In the static limit (classi- 
cal electrostatics) for the intragrain motion the Coulomb 
interaction is reduced to the effective charging energy 
Eij interaction between the total excess charges of the 
grains. Accounting for tunneling yields the screened form 
of the charging energy interaction^, which is sufffcient for 
studying the intergrain transport. We will see, however, 
that coordinate-dependent interaction modes inside each 
grain arising from the intragrain motion will be necessary 
to get a correct frequency dependence of the classical Hall 
resistance Rh of a single grain. 
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A. Kubo formula 

The conductivity in a homogeneous external electric 
field is calculated using the Kubo formula^ in Matsubara 
techniquei^: 

aab(c^) = 2e2a2-d J_ [u^^{co) - U^biO)] (3.8) 



HabH =^nab(c^,i-j), (3.9) 
j 

where 

nab(c^,i-j) = / ^ dTe^"^(r,/i,a(T)/j,b(0)) (3.10) 

is the current-current correlation function, 

/i.a(T) = Xi+aa(T) - Xu+^{t). (3.11) 

Here uj £ 27rTZ is an external bosonic Matsubara fre- 
quency (Z is the set of integers), a and b are the lattice 
unit vectors. The factor 2 in Eq. (jS.Sp stands for the spin 
degeneracy coming from one electron loop. The vector 
a denotes the direction of the current component and 
b points along the external electric field that causes the 
current. For example, if b = Gj,, then a = e^r for Hall con- 
ductivity ffxy = '^e^By and a = Gj, for longitudinal con- 
ductivity ayy = (JeyBy ■ Further, A{t) = e^'^Ae~^'^ is the 
Heisenberg operator in Matsubara technique. The oper- 
ator of the tunneling current through the contact con- 
necting the grains i and i + a actually equals —ieIi^a{T)i 
we extracted {—ie)^ from nab(w, i — j) for further conve- 
nience. The average (. . .) in Eq. p.lOp implies both the 
quantum mechanical thermodynamic averaging with re- 
spect to H and averaging over the intragrain and contact 
disorder according to Eqs. (|3.3p and (|3.6p . The contact 
between the neighboring grains i -f- a and i will be further 
identified by the pair (i -I- a, i). 

The correlation function Hab(w,i — j) represents the 
current running through the contact (i -I- a, i) in response 
to the voltage applied to the contact (j -I- b, j) only. The 
sum over j in Eq. (j3.9p means that the contributions from 
all contacts have to be considered. 



IV. TECHNIQUE 

The current-current correlation function n^bi'-^, i ~ j) 
is calculated using diagrammatic technique. Let us first 
discuss its details neglecting the Coulomb interaction He 
[Eq. (1221)] in H [Eq. ^1^] completely Technically, for a 
given pair (i -I- a, i) and (j -I- b, j) of contacts one expands 
Eq. (p^Ol) both in the disorder potential ;7(ri) [Eq. (l3?2l) ] 
and tunnelling Hamiltonian Hf [Eqs. (|3.4p and p.Sp ]. 
Each diagrammatic contribution to nab('^, i~j) is a loop 
of two Green functions connecting the contacts (i + a, i) 




FIG. 3: Different types of diagrams for the current-current 
correlation function nab(t^,i — j) [Eq. (|3.10p ] neglecting 
Coulomb interaction [Eq. (|3.7p ]. Diagrams of types (a) and 
(b) that contain oscillating at Fermi wavelength Xf = Pp^ 
functions in coordinate representation vanish after the inte- 
gration of the contacts surfaces give 0. In diagram (a) two 
different contacts are connected by a single Green function 
G{e,si,S2); in diagram (b) two different contacts are con- 
nected by two Green functions C/(e, si, S2) and 5(e, S2, si) hav- 
ing the same (sign of) energies, (c) The only type of "allowed" 
diagrams that do not contain oscillating functions and give 
nonvanishing contributions: the two contacts (i -|- a, i) and 
(j + b,j) with external tunneling vertices (wavy lines) are 
"capped" by the Green functions G(e, s, s) from one of their 
sides and connected by two Green functions Q{e + uj, Si+a, sj) 
and Q{e, sj, si+a), the "paths" of which through other contacts 
coincide. For energies, such that {e + u>)e < 0, the diffuson D 
[Eq. (|4.1|l ] in each grain along this path arises, (d) Diagram 
for the longitudinal conductivity a^J — a^~'^GT [Eq. (|1.4|l ] in 
the leading order in pt/so <C 1. 



and (j + b, j). Then one averages this loop over the in- 
tragrain and contact disorder according to Eqs. (|3.3p and 
(ISH). 

Of course, many different possibilities of drawing such 
loop can be considered (see Fig. 131). However due to the 
general properties of the Green functions in the coor- 
dinate representation and the assumption employed in 
Eq. (|3.5p that tunneling occurs from the vicinity of the 
contacts, but not from the bulk of the grain, a lot of them 
can be ruled out even before averaging over U{r). 

Consider a Matsubara Green function C/(£, r,r') of 
an arbitrary grain for a given realization of the disor- 
der potential U{t). The Green function t/(e,r,r') cx 
^ipp\r-r |sgn£ oscillatcs at the Fermi wavelength \p — 
as a function of the difference r — r'. Since we assume 
the grain size a and the size of the area of the contact 
much greater than Xp, this fact excludes the following 
possibilities. 

(i) If two different contacts Si and S2 are connected 
by a single Green function t/(e, 81,82) in a given grain 
(Fig. 13^), then integration over the contacts surfaces 
J dsids2G{e,Si,S2) gives zero due to the rapid oscilla- 
tions of the integrand. 

(ii) If two different contacts Si and S2 are connected 
by two Green functions ^(£,81,82) and f/(£',S2,Si) (or 
^/(e', Si, S2)) in a given grain having the same signs of en- 
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ergies, ee' > 0, (Fig. [8)3) then, again, their product is an 
oscillating function and J dsids2G{£,Si,S2)G{£' ,S2,Si) 
also gives zero. 

So, the only objects of the diagrammatic technique 
that "survive" inside the grains are those that do not 
contain oscillations at the Fermi wavelength Af in their 
coordinate dependence (FiglSJ;). These are: (1) the 
single Green function Q(e,s,s) with coinciding coordi- 
nates s on the contact surface; (2) the product of two 
Green functions with pairwise coinciding coordinates 
and opposite signs of energies: Si, S2)C/(£', S2, Si) or 
Si, S2)C/(£', Si, S2) with ee' < 0. After disorder av- 
eraging such products of two Green function give well- 
known electron propagators for a single isolated grain: 
the "diffuson"^ 

Diu;,r,r') = -L{g{e+Lo,r,r')g{e,r' ,r))u , ie+Lu)e < 0. 

(4.1) 

and the "Cooperon" 

C(w, r, r') = ^{g[e+uj, r, v')g{e, r, r'))c/, [e+uj)e < 0. 

Cooperons will be important for weak localization 
effects^^, while in this paper we only consider the dif- 
fuson ()4.ip in detail. 

We are left with the following general type of diagram 
(in the absence of Coulomb interaction and weak local- 
ization effects) for nab('-^, j) shown in FiglSJ:: (1) each 
contact (i-l-a, i) and (j + b, j) with external tunneling ver- 
tices must be "capped" by the Green function G{e,s,s) 
from one of its sides (one cannot "construct" a diffuson 
from G(e, s,s), since only one energy e is available, see 
Fig. [3)3); (2) two Green functions g{e + tj,Si+a, sj) and 
g{e, Sj, Si+a) connect the contacts (i -I- a, i) and (j + b, j) 
from the opposite sides and their "paths" through dif- 
ferent contacts must coincide. Therefore, in each grain 
along this path the diffuson D [Eq. (|4.ip ] of this par- 
ticular grain arises. The arising product of two Green 
functions with pairwise coinciding coordinates defines the 
diffuson of the whole granular system 

Viiv, r;, = -L{g(^s+Lj, n, rO^(e, rUi))r/.t, ie+io)e < 

(4.2) 

Contrary to Eq. (|4.ip . the points rj ri may belong to 
different distant grains i and j now. Each diagrammatic 
contribution to T> [Eq. (|4.2p ] is factorizcd into the prod- 
uct of intragrain diffusons D [Eq. (|4.1[) ] connecting differ- 
ent contacts inside the grain and tunneling probabilities 
expressed via the tunneling escape rate F. 

To get the conductivity crab(w) one should sum 
nab('^, i ~ j) over all j according to Eq. (|3.9p . An impor- 
tant observation is that due to this summation the in- 
tragrain diffusons always enter the expression for CTab('-^) 
as a difference D{lu, 81,82) — D{lj, 83, 84) of the diffusons 
connecting different contacts Therefore the zero mode 
1/{\lj\V) (see Sec. IIV Al below) drops out and the con- 
tribution to CTab(w) comes from nonzero modes with "ex- 
citation energies" of the order of the Thouless energy 



i?Th^^- Each pair "grain -I- contact" brings a factor 
F/i?Th oc gr/go, given by the ratio of the tunneling con- 
ductance gT to the conductance of the grain go- 

What does the above procedure amount to? It ap- 
pears, that this procedure reproduces exactly the solution 
of the classical electrodynamics problem for the conduc- 
tivity of a granular medium, provided each tunnel contact 
is viewed as a surface resistor with conductance Gt'^^- In 
principle, this approach allows one to study both LC and 
HC of the granular system for arbitrary ratio gr/go- For 
example, the classical formula 

^(0) ^ GtGq 



for LC can be obtained (the contact Grp and grain Gq 
resistances connected in series). Its expansion 

aj°] = o'-^Gt - G|/Go + Gl/Gl - . . .) (4.3) 

in gx / go corresponds to the expansion of the diffuson 
V in the intragrain diffusons D. Each subsequent term 
in Eq. (|4.3p corresponds to including contacts (j -I- ea;, j) 
more and more remote from {i + e^,!) in Eq. (j3.9p . 

However, for the system with well-pronounced granu- 
larity [gx <C go, Eqs. (|1.2p and (|1.3p ] one does not need to 
sum the contributions from all distant contacts (j -I- b, j) 
in Eq. 

It is sufficient to consider the lowest nonvanishing order 
in gx/go <C 1, given by the closest contacts. In fact, 
for LC a^J (a = b = e^^) considering non-zero-mode 
intragrain diffusons is not necessary at all, since the first 
term Gt [Eq. p.4p ] of the expansion (|4.3p is obtained 
from a single contact (j — i) without expanding Eq. (|3.10p 
in Ht (see Fig. [3ji). Including the closest contacts (j = 
i, i ± e^; in Eq. ()3.9p ) via the intragrain diffusons D will 
give the next term —G^/Gq in Eq. ()4.3p . which is a small 
correction to Gt- 

On the contrary, for the Hall conductivity axy (a — 
ex,h = ey) the expansion in gr/go starts from the term 
G^Rh [Eq. (HH])] analogous to -G|/Go in Eq. To 
get this term one must connect the contacts (j -I- ej,, j) 
in the y direction closest to the contact (i -|- e^:, i) in the 
X direction via the intagrain diffusons D [i.e., take into 
account the terms with j = i — Gj^, i, i + e^,, i -|- e^, — Gj, 
in Eq. (|3.9p ]. Thus, considering nonzero diffusion modes 
for Hall transport is inevitable. 

The above considerations also explain why expanding 
in the tunneling Hamitonian Ht is a "legal" procedure 
in the metallic regime, even though the dimensionless 
tunneling coupling constant (77^ ^ 1 is large. The answer 
is that the actual expansion parameter is the ratio gr/go- 

Before we proceed with the Hall conductivity, we con- 
sider important building blocks of our diagrammatic 
technique: the intragrain diffuson in the presence of mag- 
netic field and the screened Coulomb interaction. 
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e + LO 
D 



e + u 



Do 



+ Do 



e + u! 



D 



FIG. 4: Diagrammatic representation of the integral equa- 
tion (14. 4p for the diffuson D{ijj, r,r') (gray block), defined by 
Eq. (14. ip . Fermionic lines stand for the disorder-averaged 
Green function G{e, r, r'), dashed line denotes the correlation 
function ()3.3[) of the random potential. 



A. Intragrain diffuson 

The diffuson D{Lu,r,r') of a single isolated grain is de- 
fined by Eq. (|4.ip . The averaging over the disorder po- 
tential C/(r) in Eq. (14. ip is done using the conventional 
diagrammatic techniquc^^. In the "noncrossing" approx- 
imation valid for weak disorder {ppl 3> 1, where I = vpTQ 
is the electron mean free path and vp is the Fermi veloc- 
ity) the diffuson is given by the sum of ladder-type dia- 
grams. This standard series can be expressed in terms of 
the solution of the integral equation (Fig. |4]) 

D{Lu,r,r') = Do{Lu,r,r') + - [ dxDo{LJ,r,jc)D{uj,x,r'), 

To J 

(4.4) 

where 

Doito, r, r') = -^G{e + uj, r, r')G(e, r', r) 

is the "ladder step" and G{e+Lj, r, r') = {Q{e+uj, r, r'))u 
is the disorder-averaged Green function of the grain. 

In the diffusive limit [l <^ a, lutq <C 1) the integral 
equation (|4.4p can be reduced to the differential diffusion 
equation (we assume w > from now on) 

{lu - DoVl)D{Lu,r,r') ^ 6{r ~r') (4.5) 

where Dq = vpl/S is the classical diffusion coefficient in 
the grain {Dq is not affected by magnetic field, such that 
w//ro ^ 1). Equation (|4.5p must be supplied by a proper 
boundary condition at the grain surface, which we derive 
in Appendix [Xl in the presence of magnetic field: 



(n, VrD)\reS = (^HToit, VrD)\reS- 



(4.6) 



Here n is the normal unit vector pointing outside the 
grain, t = [n, H]/H is the tangent vector pointing in the 
direction opposite to the edge drift. Equation (|4.6p is 
due to the fact that the current component normal to the 
grain surface vanishes, its RHS describes the edge drift 
caused by the magnetic part of the Lorentz force. Note 
that only due to the boundary condition (|4.6D the diffu- 
son D "knows" about the magnetic field. All information 
about the magnetic field in the system is now contained 
in this boundary condition and the nonzero HC we will 
obtain is due to the nonzero RHS of Eq. (|4.6p only. The 
main consequence of Eq. (|4.6p . crucial for the Hall effect, 
is the directional asymmetry 

D{uj, r,r') ^ D{uj,r',r) 



v{n) 




FIG. 5: Diagrammatic representation of the integral equa- 
tion (|4.10p for the screened Coulomb interaction V{Q,ri,r'j) 
[wavy line, see Eqs. (|1IT1I,(|32SJ,(|1IIS1) below]. Zigzag line 
represents the bare Coulomb potential Vc(r — r') — e^/|r — 
r'l and electron loop the polarization operator 'P(r2,ri,rj) 
[Eq. (BUJ]. 



of the diffuson for H ^ 0. For = 0, Eq. (gH) reduces 
to the well-known Neumann boundary condition. 

The solution to Eqs. (|4.5p and (|4.6p can be presented 
in the form 



^(-,r,r') = -^ + E 



n>0 ' 



(4.7) 



where </>„ are the eigenfunctions of the problem 

-Vr0„ = ql<j)n, (n,Vr0„)|s = ^ff Tq (t , Vr0„ ) | S , 

7„ — -Do9n is "diffusion spectrum" , and V is the 
grain volume. The functions 0„ satisfy the orthonormal- 
ity condition 



dr0* (r)0„/(r) = 5„ 



(4.8) 



There always exists a uniform solution (j)Q{r) = 1/VV 
with the zero eigenvalue 70 = giving the zero mode 
l/{u!V) in Eq. (|4.7p . The lowest excited mode 71 ~ 
E'Yh = Dq/q"^ defines the Thouless energy scale i?Th • The 
zero mode l/(a;V) describes the fact that at time scales 
much greater than the traversal time l/-ETh the proba- 
bility density to find an electron is distributed uniformly 
over the grain volume. Information about nontrivial in- 
tragrain dynamics is contained in nonzero modes: 



^0 ^ + ^" 



(4.9) 



We will see that for Hall effect, for which the intra- 
grain dynamics is essential, only the non-zero mode part 
D{uj,r,r') of the diffuson D{uj, r,r') enters the expres- 
sions for HC and HR, whereas the zero mode 1/{ujV) 
simply drops out. 



B. Screened Coulomb interaction 

Within the random phase approximation the screened 
Coulomb interaction is given by a diagrammatic series, 
that can be obtained as a solution of the integral equation 

(Fig. ED 

V{n, r;, r^) = K(ri - r]) - ^ / ^XkdxJ x 

k,l 

xVcin - Xk) i^vin, xk, xO v{n, xi, r'.), (4.10) 
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where il S 27rTZ is a bosonic Matsubara frequency and 
Vc(r — r') = e^/|r — r'| is the bare Coulomb interaction 
Eq. (|3.7p . Just hke for ordinary disordered metals the 
polarization operator of the granular system is defined 
as an electron-hole loop 



(4.11) 

(2 comes from the spin degeneracy) and can be expressed 
in terms of the diffuson (|4.2p of the system (we assume 
n > 0): 

V{n,rur'^ - 2 [S.^Sir, ~ - r!I?(r!, r;, rp] (4.12) 

Since the Coulomb potential Vc{r — r') satisfies the Pois- 
son equation 

-V?K(r - r') = ATre^6{r - r') 
Eq. (|4.10[1 can be rewritten in a differential form 

-Y, f dxk7'(r!,ri,Xk)F(17,Xk,r(), (4.13) 

where is the Debye screening radius, l/r|, = A-kc^v. 

Depending on the approximations used for T), one ob- 
tains different forms of the screened potential V(Jl, n, rj). 



1. Coulomb interaction V(Q.,Vi,r'^) for isolated grains 

First we obtain the screened potential 1^(57, ri,rj) ne- 
glecting tunneling between the grains. In this case the 
polarization operator (|4.12p takes the form 

7'(r!,ri,rp =5ijP(l],ri,r!), 

where 

P(r2, r, r') = 2 [5(v - r') - ^Di^., r, r')] 
is the polarization operator of a single isolated grain: 

P(r!, r, r') - PnmMrWnW), Pni^) = 2-^, 

(4.14) 

[see Eq. ()4.7I) ] . Considering the limit when (i) the spatial 
scales q^^^ ^ a ^ ro are much greater than the Debye 
screening radius ro, (ii) the frequencies fl <C cr|^ are 
much smaller than the grain conductivity crf^ (trf^ oc 
Do/rj-,), we can neglect the LHS of Eq. (|4.13p altogether. 
Following the lines of Ref. [l^, we get 



(4.15) 



Here E'y = e^(C ^)ij is the charging energy matrix of 
the granular array (Cy is the capacitance matrix, see e.g. 



Ref. l20r ). The characteristic scale of E'ij is Ec = e'^/^na), 
where k is the dielectric constant of the array. The charg- 
ing energy Eij appears from the zero mode 1/{^IV) of 
the diffuson D{n, r, r'). On the contrary, the coordinate- 
dependent part of the interaction inside the grain is due 
to the nonzero diffusion modes of D{il, r, r') and equals 



n>0 



(4.16) 

For qnTD <C 1 and fl ^ crfj^ this part is completely 
screened and equal to the inverse intragrain polarization 
operator ((4ll)) . Vni^) = l/[F„(f7)]. 



2. Coulomb interaction V{0.,ri,r'j) taking tunneling mto 
account 

Now we take tunneling into account. This modifies 
the expression for the diffuson P(il,ri,rj) in Eq. (14. lip , 
which now becomes nondiagonal in the grain indices i, j. 
Let us rewrite the diffuson V in the following form: 

V{n, n, r() = Si^D{n, n, + dD{Q, ri, rp. 

The part SD{fl, ri,r'^) is responsible for tunnehng and 
vanishes, if tunneling is absent. If we leave only the zero 
intragrain modes (OD limit) in 6D{^l, ri,r'^), the diffuson 
equals 

I?(f],ri,rj) =<5ij^(17,ri,rp + il?o(r!,i,j), (4.17) 

where D{n, r, r') [Eq. (|4.9p ] is the non-zero-mode part of 
the intragrain diffuson and 

I?o(r!,i,j)-E^"''^"'^^"("'q)' 



(4.18) 



is the diffuson for the whole granular system with OD 
limit in each grain. The "kinetic term" Fq in Eq. (|4.18p 
equals 



rq-2r^(i 



where F = 2TTi'tQSo/V is the tunnehng escape rate (S'o 
is the area of the contact), f3 = x,y for d = 2 and 
/3 = x,y,z for d = 3, q G [— 7r/a, Tr/a]'^ is the quasimo- 
mentum of the granular lattice, and the sum . . . = 

/ '(2^)^' ■ • ■ denotes the integration over the first Brillouin 
zone [— 7r/a, Tr/a]'^ . 

According to Eq. (|4.17p the polarization operator 
(|4T2)) takes the form 

r;, r'.) = SiiPin, r;, r^) + ^Vo{n, i, j). 
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where 

7'o(r!,iJ) = 2[5ij - r!I?o(f^,iJ)] , Vo{n,q) = 2^y^ 

is the zero-mode polarization operator of a granular 
system". 

Accounting for tunnehng according to Eq. (|4.17p re- 
sults in the screening of the charging energy Eij in 
Eq. (|4.15p only, whereas v{U,ri,r'^) remains unchanged. 
As a result, we get for the screened Coulomb interaction 
of the granular system: 

V{n, r;, r'.) = Vin, i, j) + 5ij v{n, n, r^, (4.19) 

where 

q 

is the screened form of the zero-mode interaction^, 

£;,(q) = ^e-^"i('--'>£;i_j, 

i 

and v{^l, r, r') is given by Eq. (|4.16p . 

The form (|4.19p of the screened interaction will be suf- 
ficient for us. We will see that the nonzero interaction 
modes v{n,r,r') inside the grain will be necessary to 
get a correct classical expression for the Hall resistance 
Rh of the grain and the screened zero-mode interaction 
V{Q,i,i) will be sufficient for calculating quantum cor- 
rections to the classical result. Significant quantum cor- 
rections to HC and HR arise from the frequency range 
51 <C gxEc {grEc is the inverse i?C-time of the pair "con- 
tact -|-grain"), when F(ri, i,j) is completely screened by 
the intergrain motion: 

V. HALL CONDUCTIVITY 

Having obtained important building blocks of the 
diagrammatic technique, the intragrain diffuson and 
screened Coulomb interaction, we can proceed with our 
main goal: calculating Hall conductivity. 

We start by neglecting Coulomb interaction 
[Eq. (P?7)) ] in H [Eq. (PTTj) ] completely. As explained 
in Sec. IIVI in order to compute Hall conductivity cr^y 
(a = ea;,b = Gy) in the lowest nonvanishing order in 
gr/go ^ 1 one has to consider the contacts (j -I- e^, j) 
in the y direction closest to the contact (i -I- Ga; , i) in 
the X direction. Calculating the current through the 
contact (i + e^,!) (denoted further Sq) one has to con- 
nect the contacts 81,82,83,84 to 89, corresponding to 
j = i — Gy, i, i -I- e^,, i -I- 63; — in Eqs. p.9p and (|3.10p . 
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FIG. 6: Diagrams giving the contribution a'^y^\Lj) [Eq. (|5.8|l ] 
to the bare (without quantum eflects) Hall conductivity 
(7'^y{u}) [Eq. (|5.14p ] of the granular metal. The contacts s^, 
a = 1, 2, 3, 4, must be connected to the contact sq by the in- 
tragrain diffusons, as shown in diagrams (a), (b), (c), and (d). 
The diagrams are offset for clarity, the contact so in each dia- 
grams denotes the same contact. For each diagram four pos- 
sibilities of attaching external tunneling vertices (wavy lines) 
must be considered, as shown in Fig. [71 only one choice is 
shown here. 




FIG. 7: For each diagram in Fig. [6] four possibilities [two for 
each contact according to Eq. (|3.1ip ] of attaching external 
tunneling vertices (wavy lines) must be considered. 



respectively, by the diffusons D(w,8o,8a), a — 1,2,3,4, 
of a single grain, as shown in Fig. [6l 

Let us consider the contribution to the cor- 

relation function Ha;y [w, i — (i — e^)] [Eq. p.lOp ] from the 
contact si. From now on we assume w > 0, the arrow 
subscript denotes the direction of the diffuson accord- 
ing to Fig. [6l the superscript "0" stands for the "bare 
value" without quantum effects of Coulomb interaction, 
the superscript "1" is introduced, since there will be an- 
other contribution "2" to HC, see Sec. IV Al According 
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to the diagram in Fig. [6l we get: 



where 



X [Gi+e, (e + uj, So, So) - Gi+e, (e, sq, So)] 
X [Gi-ej, (e + si, si) - Gi-ey (e, si, si)]. 



Each end of the diffuson Dj of the grain i is "capped" 
by the Green functions Gi+e^ and Gi_ej, of the adjacent 
grains i + e^, and i — ej^. We do not write the grain 
subscripts further. The difference G{e + to) — G{e) for 
each contact arises due to two possibihties of choosing 
the external tunnehng vertex in h a: Xi+a,i or — Xij+a, 
see Eq. p. lip and Fig. [71 For the Green functions at 
coinciding points one can use their bulk expression (with 
H^O) G-i(£,p) =*£-^(p)+ 2^sgne: 



fn,a ^ 'g2 j cJsodSa^f)„(so)(/i*(Sa) (5.6) 

with a = 1,2,3,4 for a — ,\, respectively 
(Fig. [6]). The factor /„ takes care about the geometry and 
gives a convenient compact form of the contributions. It 
will be especially helpful for studying interaction correc- 
tions to HC. We can rewrite Eq. (|5.3p as 



d^[G(e + c.,p)-G(e,p)] = 



G(£ + u, So, So) - G(e, so, so) = 

— 27rw, (e + iL>)e < 0, 
0, {e + uj)e>0. 



Therefore, we get 

n^'^^(w) = w — [ dsodsiD{ui,So,Si), (5.1) 

where gr = 27r(i^io)^'S'o is the conductance of a tunnel 
contact, 5*0 is the area of the contact, and uj arises as 

27rrE_^<,<ol 

Carrying out the same procedure for the remaining 

contacts S2,S3,S4 and paying special attention to the 

signs of the contributions, for the total contribution 

nir^-) = n5^^)(.) + n(^'^)(c.) + n^^^)(c.) + nl^'^'(c.) 

(5.2) 

to Ilxy{u!) [Eq. (|3.9p ] from the diagrams in Fig. [51 we 
obtain 

n(°/)(.;) = u;^ [DA^) - ^n(^) + ^/(^) ~ D^^)] . 

(5.3) 

Here 

Da{uj) ^ I dSodSaD{uj,SQ,Sa) (5.4) 

with a ~ 1,2,3,4 for a =/',\,j/,\, respectively 
(Fig.li. 

Using the expansion (|4.7p for the diffuson, we see that 
due to the sign structure of Eq. (|5.3p the zero mode 
l/{u!V) drops out of it. Therefore, retaining only the 
zero mode in Eq. (|4.7p would give just in Eq. (|5.3p . and 
we are forced to take all nonzero modes into account. 

According to the structure of Eq. (|5.3p we introduce 
the following auxiliary quantity: 



n(°/)(c.) = c. 



5t \ ^ fn 



V ^ UJ + 7 „ 

n>0 



(5.7) 



and the contribution to Hall conductivity (|3.8p corre- 
sponding to HS/^''(a;) takes the form 



-d.9T \ ' /" 



a(°/)(.)=2eV-'^^5: 



n>0 ' 



A. Correct tj-dependence 



(5.8) 



(5.5) 



As we show further in Sec. IV Bl the expression (|5.8p 

for afy'^{u}) at zero frequency w = reproduces exactly 
the result (|1.5p for HC of a granular medium obtained 
from the solution of the classical electrodynamics prob- 
lem. Therefore, it would be natural to expect such cor- 
respondence with classics for all lo. 

However, the obtained result (|5.8p does not agree with 
the classical formula (|1.5p at finite frequency w > 0. 
Indeed, in classical electrodynamics the resistance of a 
metallic sample itself is frequency-independent up to 
very high frequencies w ^ of the order of the grain 
conductivity crflj. (crflj. cx g^Ec is the inverse _RG-time 
of the grain)'^^. So the HC we are looking for should 

be given by the zero- frequency expression a^xy^^ {uj — 0) 
for all w ^ af^. We see, however, from Eq. (|5.8p that 

'Jx'y^\'^) has a diffusion-like dispersion at Thouless en- 
ergy (when LO ~ i?Th) which contradicts the classical pic- 
ture. 

What is yet missing in our approach? Apparently, 
one must consider the Coulomb interaction inside the 
grain. Indeed, the diffuson D{ijj,v,v')^ appearing in 
Eqs. (|5.3p . (|5.7p . (|5.8p . describes the propagation of elec- 
tron density, but it does not take into account that elec- 
trons are charged and can interact. The correct way to 
include the classical effect of Coulomb interaction is to 
insert the interaction line into the diffuson as shown in 
Fig. [SId. In its turn, this interaction is screened by elec- 
tron motion, and this can be accounted for by inserting 
polarization operators into the interaction lines, as shown 
in Fig. [5};. The summation of the resulting series in Fig. [51 
yields the screened form of the Coulomb interaction in a 
grain F(a;, r, r') = ^(w, i, i) -I- w(w, r, r') [Eq. (|¥Tg)) ]. and 

we obtain an additional to T\^'^'^\uj) [Eq. (|5.ip ] contribu- 
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FIG. 8: Complete set of diagrams for the bare (without quantum effects) Hah conductivity aij (to) [Eqs. (|5.14|l ] of the granular 
system, (a) One starts by connecting the relevant contacts by the intragrain diffusons (see Fig. |6]). (b) To get a correct 
oi-dependence one should take the Coulomb interaction into account by inserting the interaction line into the diffuson. (c) 
Coulomb interaction, in its turn, is screened by the intragrain motion and one should consider insertions of the polarization 
bubbles into the interaction line, (d) The summation of the resulting series yields an additional to Il^y^'' (u) [(a), Eq. (|5.7|l ] 
contribution nSj^-'(cij) [(d), Eq. (|5.10|l ] to the current-current correlation function. The sum 
[Eq. (|5.11|l ] gives a correct classical expression (|5.14|l for the Hall conductivity axy{i-o). 



tion 

n'^y\uj) = 2^2-^ J dsodsij dvdv' X 

So, r)F(cj, r, v')ujD{uj, r', Si). 

Here the integration with respect to r and r' is done over 
the grain volume, the spin degeneracy factor 2 comes 
from an additional electron loop. Due to the orthogonal- 
ity of the eigenfunctions (pn [Eq. (|4.8p ] the zero modes of 
D and V drop out of the needed combination 

4°/' H = n^y' (c.) + u^^ (.) + u^^ [u] + nl^'^) {.) , 

(5.9) 

and we get an additional to niy^''(w) [Eq. (|5.7p ] contri- 
bution (Fig. Hi): 



n(<;^^^)(c.) = 2^5]/„^.„(..)- 



(5.10) 



Summing the contributions (|5.7p and (|5.10p we obtain 

n(°)(c.) ^ ni°/)(^) + ni';,-^)(^) = |:^c./„r„, (5.11) 

ri>0 

where 

1 / 9,,, \ 1 

(5.12) 



1 / / X 2w 

1 + Vn{u}) 



for oj < af^ [see Eq. (|iT^ ]. Equations and ([?T^ 

lead to the final classical expression for the Hall conduc- 
tivity [Eq. dSilD]: 

<)(c.) = 2eV"'^^^:^, (5.13) 
or, going back to diffusion propagators, 
ailH^) = 2e'a^"'^{D^ ^D^+Dy- 5^), (5.14) 

where Da — f dsodSaD(so,Sa), with a = 1,2,3,4 for 
a \ , »/, \, respectively, and 

^ ^^0„(r)C(rO 



is the diffuson without the zero mode at = satisfying 
Eqs. and ([iTB]) with = 0: 



-i?oV^^(r,r') =<5(r-r'), 

(n,VrL')|s - WHT-o(t,Vr5)[s 



(5.15) 



It follows from Eqs. (|5.15p that Z3(r,r') is a Green func- 
tion for the Poisson equation. Actually the propaga- 
tor D{r,r') should not be termed "diffuson" anymore, 
since it describes the propagation of electron density with 
Coulomb interaction taken into account, i.e. the real con- 
duction process. 

Equation (|5.14|) constitutes our main result for HC 
in the absence of quantum effects. We stress that the 
diagrammatic series in Fig. [5] leading to Eq. (|5.14p de- 
scribes the classical effect: propagation of electron den- 
sity in a disordered metallic sample. The temperature- 
independent result (|5.14p is valid for arbitrary tempera- 
ture T and arbitrary size a of the grains (not necessarily 
small grains and T <C i?Th)- Temperature will be rele- 
vant for quantum effects of Coulomb interaction, which 
we consider in Sec. IVII 



Properties of Eq. \5. 14\ l for Hall conductivity a^xy ■ 



Let us discuss the basic properties of Eq. (|5.14p . For 
simplicity, we assume that grains have reflectional sym- 
metry in all three dimensions. Then — Dy and 
~ D-K^ due to this symmetry (for H ^ 0, too). 
At zero magnetic field {H — 0) we have Dy — D-^ 
and Dy — D'\^ due to the time-reversal symmetry 

D{uj, r,r') = D{ijj,r' ,r), and therefore a'^xyi'-ij,!! = 0) = 
0. The nonzero differences D 



D 



D. 



arise only due to nonzero RHS of the boundary condi- 
tion Eq. (|4.6p . which represents the edge drift. To un- 
derstand the sign of D y — and Dy — we recall 
that the diffuson D(cL',r,r') describes the probability of 
getting from point r' to point r. In nonzero field (H ^ 0) 
the edge trajectories for Dy = Dy are shorter (if e > 
is assumed) than those for = D^, and therefore 
D^-D\ =Dy-D^ >0. 
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1 1 

linear in i^hTo: one can estimate 



Since Da oc 



and the difference — is wfien tfie current ly = I passes tlirough the grain in the 

y direction. The current density 



D 



where 



Hxy 



1 



^hTq 2 P'xy 
oc e V 



(5.16) 



H 



nec 



is the specific HR of the grain material expressed in terms 
of the carrier density n in the grains [Einstein relation 
(t|^ = 2e'^iyDo was used in Eq. (|5.16p ]. We see that 
D y — does not depend on the intragrain disorder, 
described by the scattering time tq. The proportionality 
coefficient in Eq. (|5.16p is determined by the shape of the 
grains only. Thus, for HC [Eq. (|5.14|) ] of the system we 
get 



and for HR [see also Eq. pTI)) ] 



JO) 
rxy 



>xy 



(o)^2 

XX } 



oc a 



'3 

Hxy 



H 



We come to an important conclusion. The Hall resis- 
tivity of the granular system is independent of the intra- 
grain disorder and tunneling conductance. It is expressed 
solely via the carrier density n of the grain material up 
to a numerical coefficient determined by the shape of the 
grains and the type of granular lattice. 



B. Classical picture 



Let us now prove that Eq. (|5.14p for the Hall con- 
ductivity indeed reproduces the solution of the classical 
electrodynamics problem, provided one treats the tunnel 
contact as a surface resistor with the conductance Gt- 

The classical HC of the granular medium in the limit 
Gt ^ Go [Eq. (|1.3p ] can be easily presented in the form 
of Eq. (fTT5)) (see Fig. [1]). The current ly = GxVy run- 
ning through the grain in the y direction causes the Hall 
voltage drop Vh = Rnly between its opposite banks in 
the X direction, where Rh is the Hall resistance of the 
grain. Since for calculating Oxy the total voltage drop 
per lattice period in the x direction is assumed zero, the 
same voltage Vh (but with the opposite sign) is applied 
to the contacts in the x direction. Thus, the Hall cur- 
rent equals Ix — GtVh = G^RnVy which leads to the 
expression (jl.Sp for HC. 

The Hall resistance Rh of the grain is defined via the 
difference (Hall voltage) of the electric potential (/^(r) be- 
tween the opposite banks of the grain in the x direction, 



Vh = </3(Sr) - <^(s/) = Rnly, 



(5.17) 



j(r) = -o-oVr(p(r) 



XX 


x:y 





'-'xy 


^ XX 











^ XX 



(5.18) 

((To is the conductivity tensor) satisfies the continuity 
equation 



divj = q{r) 
and the boundary condition 

(n,j)|5 = 0. 



(5.19) 



(5.20) 



The charge source function g(r) is nonzero on the con- 
tacts surface only, / dsdq{sd) = I corresponding to the 
current / flowing into the grain through the contact 
and / (is„g(s„) — —I corresponding to the current flow- 
ing out of the grain through the contact s„. The sta- 
tionary form of Eq. (|5.19p is valid up to the frequencies 
Lu ~ crflj., even if / = I(t) is time-dependent, compare 
with discussion of Eq. ((1715)) in Sec. IIVBI 

Inserting Eq. ([5T^ into Eqs. ([ST^l) and ([00)) . wc find 
that Lp(r) is a solution of the following boundary value 
problem: 

-V2^ = g(r)/(7f,, (n,V^)|5-™(t,V^)|5. (5.21) 

Comparing Eq. ()5.2ip with Eqs. ()5.15p . we see that 
D{r, r') is a Green function for the problem (I5.2ip . Thus 
the solution to Eq. (I5.2ip can be written as 



ip{r) 



1 / 



dsdD{r,Sd)- / (is„D(r,s„) 



(Einstein relation crf^ — 2e^vDi^ was used). Inserting 
if{Y) in such form into Eq. ()5.17p . we obtain for the Hall 
resistance of the grain: 



Rh — 



1 



2e'^v 



D 



(5.22) 



Comparing Eq. ()5.14p with Eqs. ()1.5p and ()5.22p we see 
that Eq. (|5.14p indeed reproduces the classical result. 

This establishes the correspondence between our di- 
agrammatic approach of considering nonzero diffusion 
modes and the solution of the classical electrodynamics 
problem for the granular system. 

Luckily, for simple geometries of the grain (cubic, 
spherical) the Hall resistance Rh can be obtained 
from symmetry arguments without solving the problem 
Eq. ()5.2ip . Suppose the grain has reflectional symme- 
try in all three dimensions. Then it is clear that (1) the 
largest cross section of the grain lies in the plane of reflec- 
tion, (2) the current density j(r) is perpendicular to the 
plain of reflection at each point r of the cross section, (3) 
the absolute value of j(r) is constant on the cross section 
and therefore equal to |j(r)| = I/S, where S is the area of 
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the cross section. So, the HaU voltage Eq. (|5.17p equals 
Vh = p|y|j(r)|a — al/S. Therefore, the Hall resistance 
is 

Rh = p%a/S (5.23) 

and the Hall resistivity of the granular medium can be 
expressed in the form 

(4°j)2 n*ec 

where 

n* = a^-'^An, A = S/a^ < 1. 

The quantity n* defines the effective carrier density of the 
granular system. For a 3D sample (many grain monolay- 
ers), n* — An differs from the actual carrier density n 
of the grain material only by a numerical factor A deter- 
mined by the shape of the grains and type of the granu- 
lar lattice. For a 2D sample n* = aAn for a single grain 
monolayer or n* = dzAn in case of several monolayers, 
where dz is the thickness of the sample {dz/a is the num- 
ber of monolayers). 

We remind the reader, that Eq. (|5.24p was obtain un- 
der the following assumptions: (a) diffusive limit inside 
the grains, I <ti a; (b) the mean free path I is the same for 
all grains; (c) the tunneling conductance Gt is the same 
for all contacts. Having established the correspondence 
between our diagrammatic approach and the classical so- 
lution of the problem, we can now show that the result 
(|5.24p is actually valid in a much more general case, when 
(i) the intragrain disorder is ballistic, I > o; (ii) the mean 
free path li varies from grain to grain (iii) the tunnel- 
ing conductance GTi+a,i varies from contact to contact. 
The statement (i) follows from the fact that the above 
classical consideration leading to Eq. (|5.23p involve only 
the symmetry properties of the current distribution j(r) 
and therefore hold for the grains with ballistic disorder 
{I > a) as well. The statement (ii) is true, since the Hall 
resistance Rh [Eq. (|5.23p ] of the grain is independent of 
the mean free path li, and, hence, is the same for all 
grains i provided only their shape is the same. Finally, 
using the standard Ohm and Kirchhoff laws, we obtain 
that the Hall voltage between the opposite banks of the 
sample in the x direction depends on the total Ohmic 
current flowing through the sample in the y direction. 
Hence, the Hall voltage is essentially independent of the 
distribution of the tunneling conductances Gxi+a.,!, and 
the HR of the system is still given by Eq. (|5.24p . which 
proves the statement (iii). 

Therefore, the result (|5.24p for Hall resistivity is ap- 
plicable to real granular systems, in which fluctuations 
of the intragrain mean free path I and, most impor- 
tantly, the intergrain tunneling conductance Gt are al- 
ways present. 

To conclude this section, we obtain that the classical 
Hall resistivity [Eq. (|5.24p ] of a granular system in the 



metallic regime, being independent of parameters that 
describe Ohmic dissipation, possesses a great deal of uni- 
versality, reminiscent of that in ordinary disordered met- 
als [Eq. (fLTl) ]. 

Being classical, however, Eq. (|5.24p describes the be- 
havior of the Hall resistivity at high enough tempera- 
tures, when quantum effects can be neglected. At suf- 
ficiently low temperatures quantum effects of Coulomb 
interaction and weak localization set in and can signifi- 
cantly affect electron transport. In the next section we 
study quantum corrections to the obtained results (|5.14p 
and (|5.24p due to the Coulomb interaction between the 
electrons. 



VI. QUANTUM EFFECTS OF COULOMB 
INTERACTION 

Diagrammatic approach enables us to incorporate 
quantum effects of Coulomb interaction on the Hall con- 
ductivity into the developed scheme. We perform calcu- 
lations to the first order in the screened Coulomb inter- 
action with the expansion parameter l/gx- We assume 
the diffusive limit for the intragrain dynamics and neglect 
the diagrams that are small in ToETh oc {l/a)"^ <C 1. The 
ballistic limit can be treated similarly, although in this 
case one has to take such diagrams into account. 

Technically, one considers the diagrams for the "bare" 
conductivity cr^J shown in Fig. [8] and connects differ- 
ent electron lines by the interaction lines correspond- 
ing to the screened Coulomb interaction V{fl,ri,ry) 
[Eq. (|4.19p ]. It is important that for the quantum interac- 
tion corrections the zero-mode part V{^1, i, j) [Eq. (|4.20p ] 
of the interaction V{fl, ri, Tj) does not drop out and gives 
a contribution larger than the nonzero intragrain modes 
Sijv{n, Ti, Tj) [Eq. (|4.16p ] (we provide an estimate below). 
Therefore for the interaction lines that describe the quan- 
tum corrections to the classical result we can use the 
zero-mode part V{fl,i,j) of the interaction. Further, de- 
pending on the sign structure of energies of the Green 
functions involved, some interaction vertices are renor- 
malized by the diffusons and some are not. 

Two types of diagrams can be identified: (i) the in- 
teraction V{n,i,j) connects different electron loops of 
the diagrams like in Figs. [HI^a),(b); (ii) the interaction 
V{fl, i, j) connects points on the same electron loop, like 
in Figs. [9l(c),(d). It is straightforward to show that the 
former possibility (i) always gives zero: in each case con- 
tributing diagrams cancel each other identically , an ex- 
ample is shown in Fig. [9Ka),(b). So, we come to an im- 
portant simplification: electron loops in Fig.[8]are renor- 
malized by the interaction independently. 

The temperature T, being irrelevant for the single- 
particle classical transport [Eq. (|1.5p for a^xy], becomes 
important for quantum effects of Coulomb interaction. 
An important energy scale here is the tunneling escape 
rate F. For T > F the thermal length = ■\JD*/T < a 
for the intergrain motion {Dq = Va? is the effective diffu- 
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FIG. 9: Quantum corrections to Hall conductivity of a granu- 
lar metal due to Coulomb interaction. Important point: elec- 
tron loops in diagrams of Fig. |S] are renormalized indepen- 
dently: the diagrams with interaction line connecting differ- 
ent electron loops cancel each other, as is the case, e.g., for 
the diagrams (a) and (b). Nonvanishing contributions come 
from the independent renormalization of electron loops by the 
interaction, like in diagrams (c) and (d). 



sion coefHcient) does not exceed the grain size and only 
the contributions coming from spatial scales of the order 
of the grain size a are significant. At T < F the contri- 
butions from the scales > a exceeding the grain size 
also become important. 

We start with the former regime T ^ F in the fol- 
lowing subsection. Note that the corrections that will be 
discussed in Sec. IVI Al are specific to granular systems 
and absent in HDMs. At the same time they, as we find, 
govern the T-dependence of HC and HR in a wide range 
of temperatures. The corrections analogous to those in 
HDMs ( "Altshuler-Aronov" corrections) arise from large 
spatial scales a), are relevant at T ^ F and will be 
addressed afterwards in Sec. IVI Bl where we consider the 
case T < F. 



A. "High" temperatures T > F 

First, we consider the range of temperatures T ^ F 
greater than the escape rate F. In this regime each tun- 
neling process brings a small factor F/T ^ 1. Therefore 
the main contribution comes from the diagrams which 
contain minimal number of hops between the grains as 
compared to the diagrams for the bare HC cr^J . This 
means that the screened zero-mode interaction can be 
taken in the form [sec Eq. (|4.20p ]: 



F(f],q) 



Eciq) 



l + 2E,iq)Tj{S\n\) 



(6.1) 



and for the diffusons renormalizing the interaction vertex 
we can neglect tunneling completely, i.e., take them as 



It is very important that, since the interaction V(ri, i,j) 
is coordinate-independent within each grain, the intra- 
grain diffusons renormalizing the interaction vertex con- 
tain only the zero mode l/(|r2|V), whereas the nonzero 
modes drop out automatically due to the orthogonality 
condition ()4.8|) (we do not simply neglect them): 

J dr[Din,n,v'^vin,i,j) = ^^vin,i,i). 

Assuming T 3> F, we do not assume the temperature 
T or the frequencies much smaller than the Thouless 
energy i?Th in this section. As we will see, the only way 
to clearly identify the physics of the contributions and 
obtain a correct upper cut-off for the logarithmically di- 
vergent quantities is to include the range T, i7 ^ E'xh into 
consideration. 

As explained above, we may renormalize electron loops 
shown in Fig. [8] independently of each other. There are 
three different types of electron loops in Fig. \8\ 1) the 
(tunneling current)- (tunneling current) correlator (//) of 
the diagram in Fig. UKa); 2) the (tunneling current)- 
density correlators (In) in Fig. [5]Jb) connected by the 
screened interaction v{uj,r,r'): 3) the density-density 
correlators (nn) [i.e. the intragrain polarization opera- 
tor P(w,r,r'), Eq. (|4l4)) ]. in Fig. ^c), which are the 
insertions into the bare Coulomb interaction line. 

Since the geometry factor /„ [Eqs. (|5.5p and ()5.6|1 ] is 
the same for all diagrams and the properties of the granu- 
lar array are assumed the same in the x and y directions, 
let us draw the diagrams in the "longitudinal geometry" 
(see Figs. ll0ll2ll3p . For our purpose it is only important 
now that there is a "central" grain with a non-zero-mode 
diffuson and there are two "adjacent" grains. External 
tunneling vertices are attached to the contacts between 
the central and adjacent grains. 

For each diagram one has to take 4 possibilities (2 for 
each contact) of attaching external tunneling vertices into 
account, as in Fig.[Sl Only one such possibility is shown 
in Figs. I10I12I13] Further, for each diagram one has to 
consider (i) the up/down reversal, if the diagram does not 
transfer to itself, (ii) the left/right reversal, if the diagram 
does not transfer to itself. We introduce "up/down rever- 
sal" and "left/right reversal" multiplication factors a"/'^ 
and a'/"^ correspondingly: a"/'^,Q;'/'^ = 1, if the reversal 
is not possible, if the reversal is possible. 

The summation region over the fermionic frequency e 
of the electron loop and bosonic frequency fl carried by 
the interaction line is determined by the analytical prop- 
erties of the Green functions involved. After the integra- 
tion over the Green functions momenta the expressions 
become independent of s and the summation over s can 
be performed. This always results in the sum 

27rT F{uj,fl) + 2nT ^ F{uj,n) ^ 



0<f2<w. 

-n<e<6 



-a)<£<0 



= J2 r!F(cj, r!) + ^ F(cj, r!) ^ ^ 0^^oF(c^, i|^.2) 
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standard for the first-order interaction corrections 
calculations^. We have introduced the function 

for compactness. The diagrams considered here may 
contain either one or two summations (|6.2p . We intro- 
duce the "sum" multipHcation factor correspondingly: 
a'' = 1 or a*^ = 2. So, each "topologically unique" dia- 
gram comes with an overall multiplication factor 



i 


Expression for Xn\uj,^l) 






1 /v 

a/ 




. + + 

1 11 

VoiQ) 

uj + -ynUJ + il + 'y„u 
1 1 

{lu + 'y„y^ uj + n + 'y„ °^ 


\ 


2 


\ 


d 


1 


2 


1 


e 


1 


2 


2 


f 


1 


2 


1 



TABLE I: Corrections to the intragrain polarization operator 
Pn = 2[1 - lj/{lu + 7„)] [Eq. gH]. 



a = a'' a"/dQ,i/'-. (6.3) 

We start by considering the corrections to the intra- 
grain polarization operator. 



1. Corrections to the intragrain polarization operator 
P{lo, r,r') 

The set of diagrams for the first-order interaction cor- 
rections to the intragrain polarization shown Fig. [10] is 
the same as the one for a bulk system^ The crossed re- 
gion in diagram (c) is a Hikami bojs^ shown in Fig. [TTl 
which analytical expression for the case of coordinate- 
independent interaction is —(uj + Q + jn)- 

We present the correction to the n-th mode P„(a;) = 
2[1 — uj/{uj + 7„)] of the polarization operator P{u!, r, r') 
[Eq. (|4.14|) ] in the form 

SPi'\cu) = 2T^ 0^,na.AW(w,f)), i = c,d,e,f, (6.4) 

where the expressions for Xn\uj,fl) and the multiplica- 
tion factors ai = a^^'^a/'^ [Eq. (|6.3p ] are given in Ta- 
ble U (the factor 2 stands for spin degeneracy, the dia- 
grams are labelled in correspondence with the diagrams 
for SIl^xy'^\uj) and SIl'^y'^\Lu) below). In Table U and 
Fig. [TOl the quantity Vb(ri) = V{n, i, i) is the zero-mode 
interaction in a given grain. Summing the contributions 
(16. 4p . we obtain that the total correction to the intra- 
grain polarization operator P„(a-') due to the zero-mode 
interaction Vo(ri) vanishes: 

SP^iico) + 5P^{u) + 5P:{u) + SP'Jlo) = (6.5) 

This is an expected result, since due to the gauge in- 
variance the constant interaction potential cannot affect 
physical quantitie o^^'^^ , expressed in this case by the 
density-density correlation function. 

As a result, we obtain that the screened non-zero- mode 
Coulomb interaction ?;(w,r, r') [Eq. (|4.16p . double wavy 
line in Fig.[5Ji] does noi acquire any correction. Therefore 
we should only renormalize the electron loops (//) and 
{In) shown in Figs.[S^,d explicitly. 



2. Interaction corrections to niy''''(a;) and Iliy'^'' (a;) . 

Now we renormalize the electron loop (//) of nS/^''(a;) 
[Eq. (1121), Fig. Hi)] and two loops (/n) of nl''/^(w) 
[Eq. (|5.10p . Fig. ISJi]. The nonzero modes v„{u!) of 
the screened intragrain interaction (double wavy line in 
Fig. ISji) are not renormalized according to the result of 
the previous subsection. 

All corrections to lliy'"'^'' (a>) and ni°'^-'((jj) may be pre- 
sented in the form: 

<5n«(^) = ^ E ^ E ^--^ (6.6) 
n>o n>a 

The geometry factor /„ [Eqs. (|5.5p and (|5.6p ] arises, when 
corrections to all diagrams in Fig. l6l from the four closest 
contacts are taken into account. The sets of diagrams 
giving corrections to Ti^xy^^{oj) and nS/^''([j) are shown 
in Figs. [121 ^^nd I13i and the corresponding expressions 

for \n\ijJ,^) and ai — a^a^^'^a/'^ are given in Tables HIl 
and mil respectively. For each diagram in Figs. [T^ and 
ITSlone must take four possibilities [two for each contact 
according to Eq. (|3.1ip and as shown in Fig. IS] of at- 
taching external tunneling vertices into account. In the 
expressions, 

Fo(i]) -F(f^,i,i) = E^(^^.q) 

q 

is the "on-cite" interaction of the grain, 

Vi{n) = V{n, i + e„ i) = E cos{q,a)V{n, q), i = x,y, 

q 

is the interaction between the neighboring grains, and 

V2{ft) = V{n, i+ex+e.y, i) = E cos(gxa) cos{qya)V{il, q) 

q 

is the interaction between the next-to-the-nearest- 
neighboring grains. The interaction V{il,q) is given by 

Eq. dnn). 

Note that the diagram (d) in Fig. [T^l giving a correc- 
tion to Il^xy^\uj), does not have an analog for 11^/'^-' (w), 
because this would mean connecting different electron 
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FIG. 10: Diagrams for the interaction corrections to the intragrain polarization operator Pn{u}) = 2[1 — 01/(0; + 7„)] [Eq. (|4.14|) ]. 
The corresponding expressions are given in Table HI Gray blocks depict nonzero diffusion modes, rendered with lines blocks 
depict zero- mode diffusons 1/Q renormalizing interaction vertices, dashed lines stand for the impurity correlation function (|3.3|l . 
The crossed block in diagram (c) is the Hikami box (see Fig. Ill|) . 



-(w + + 7,0 





+ 





' I 
I 

— *- 



FIG. 11: Hikami box. Analytical expression for the Hikami 
box in case of coordinate-independent interaction potential is 

— (0 + + 7n). 
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Expression for X^\uj,Q,) 


Qf 


u/d 


l/r 


al 


— ^^vom 

Lu + -y„ il'' 

LU + 7,1 il^ 


2 


1 


2 


bl 


1 


2 


2 


cl 


(- + " + -^^")(o + 7.)^n^^°(") 

^ + n + j„n^^'^^'^ 

^ ^ Vi(n) 

+ 7,1 + n + 7„, n 


1 


2 


1 


d 


2 


1 


1 


el 


1 


2 


2 


fl 


1 1 

(o + 7„)2a; + n + 7„''^'^"^ 


1 


2 


1 



TABLE H: Corrections to (w), [Eq. ((53, Fig. [Ha)] 



loops by the interaction line, which gives 0, as discussed 
before. The crossed region in diagrams (cl) and (c2) is 
the Hikami box (Fig. [TT|) . 

Let us perform partial summation of the contributions 
(16. 6p as follows: 

a^iX^^LO, n) + a^2Xniio, O) = -4^^ ^o(0)r„, (6.7a) 



abiA^^(w, n) + ab2A^"(a;, O) = A-^Vi{n)rn, (6.7b) 



Expression for Ai'-'(o,r2) 



Lj + 7„ Sl^ 
1 1 



2o; 



o + 7„ 
a; + ri + 7 



Vi(n) 



o + 7„ 
2aj 



V^.(J7) 



o + 7„ 
2uj 



(o+7„)2n2 

1 11 

a; + 7nOi + Sl + 7„r2 
1 1 



2UJ 



LJ + n + JnO.^ 
1 1 



O + 7r 
2llI 

Vn 



-Vo{n) 



{UJ + 7„)2 LJ + n + 7, 

— 3^ ^ -Voifl) 

o; + 7n 01 + ^2 + 7n 12 



2iu 



o + 7„ 
2uj 



O + 7n 



u/d l/r 



TABLE HI: Corrections to Hi°/'(a;), [Eq. (tCTJjl . Fig. Eld)]. 
5] a.A«(c.,0) - 2 ^ J-Vb(0)-4-^Fo(0)r„, 



i=cl,c2, 
fl,f2,f3 



(6.7c) 



i— el, 
e2,e3 



.,A«(u;,0) 



1 1 

W + + 7„ 02 



P2(0), (6.7d) 



a; + O + 7„ O^ SZ^ 



(6.7e) 

where ?'„ — l/7„ is the "resistance" mode arising as a 
sum of the series shown in Fig. [5] and given by Eq. (|5.12p . 

We see that two functionally different forms arise. One 
contains the frequency-independent resistance modes 
r„ = l/7n, the other contains nonzero diffuson modes 
1/{lu + + 7„). Summing these two independently, 
we present the total first-order interaction correction 
6Ilxy{(^) to the classical result ()5.1ip for the current- 
current correlation function Ill^y{uj) as: 



(6.8) 
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^ + 7n ^ + ^ + "/n ^ + 7n UJ + Q + Jrt. 1^ + 7^ ^ + 7n. CJ + O + 7„ U + 7^, 



FIG. 12: Diagrams for the Coulomb interaction corrections to Iliy (uj) [Eq. (|5.7p . Fig. |8ja)] describing renormalization of 
the (tunneling current)-(tunnehng current) correlator {//) of Fig. (SJ^a). Open circles denote tunneling vertices placed at the 
contacts. Other elements are explained in caption to Fig. 1101 




FIG. 13: Diagrams for the Coulomb interaction corrections to niy'^-'(Lj) [Eq. (|5.10|l . Fig. |8jd)], describing renormalization of 
the (tunneling current)-density correlators (In) of Fig. [8jd). 



where 



is obtained from Eq. 



according to Eq. 



(6.9) 



1 



n>0 n>o 



and 



n>0 O>0 ' 

x-^[V2{n) + Voifi)-2Viin)]. (6.10) 



The corresponding correction Saxy{uj) — 5a'^^{uj) 
^'^xv'i^) to Hall conductivity 



3.11) 



.12) 



Let us discuss the obtained results (|6.8p . (|6.9p . (|6.10p . 
First of all, note that the corrections 6H'^^{uj) and 
SU^^{iLj) cannot be attributed to certain separate sets 
of diagrams: the diagrams (al),(a2),(bl),(b2) contribute 
to SU'^^iuj) only [Eqs. ([67a|) and ((67H)) ]. the diagram 
(d) contributes to SU^^{lj) only [Eq. (j6.7dp ]. but the 
rest of the diagrams contain parts corresponding both 
to dll'^^{uj) and (511^/^ (w), as seen from Eqs. (|6.7cp and 

Note also that the corrections Sli^^{uj) and 611^^ (uj) 
vanish separately in the case of constant interaction po- 
tential, i.e., when Vo{fl) = Vi{n) = V2{^). This property 
is not accidental and is enforced by the gauge invariance 
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(see, e.g., Refs. [2lll23 ): constant interaction potential re- 
sults in a shift of the chemical potential of the whole 
electron system and therefore does not affect physical 
quantities expressed diagrammatically as chains of closed 
electron loops (see Fig. |8]). This fact serves as an im- 
portant check of our resufts ^M,^M, and KTO\i . We 
remind the reader that Eq. (j6.5p is also a consequence of 
the gauge invariance. 

The physical processes leading to the corrections 
SH'^^{uj) and SU^^{lu) are most clearly identified from 
the diagrams that contribute to either one of the quan- 
tities, i.e., from the diagrams (al),(a2),(bl),(b2) for 
SH'^^{uj) and from the diagram (d) for SU]^^{lj). 

3. "Tunneling anomaly" contribution Sli^^yis^j) 

First, consider the correction 5lC^^{bj) [Eq. (|6.9p ]. 
Clearly, the diagrams (al),(a2),(bl),(b2), contributing 
solely to 5\i^y{ui), describe the effect of Coulomb inter- 
action on the process of electron tunneling through the 
contact. Therefore, the correction 5Ii^y{ijj) should be 
attributed to the tunneling anomal'i^^'^ (TA) effect in 
a granular metal. The correction SIl^^{uj) corresponds 
to the independent renormalization of the tunneling con- 
ductances Gt of the individual contacts in formula ()1.5|) 
for the bare classical HC a^J oc G^Rh, whereas the Hall 
resistance Rh of the grain remains unaffected. 

Indeed, the r2-independent resistance modes r„ 
[Eq. (|5.12p ] can be taken out of the sum over in 
Eq. (|6.9p . After that the summed over expression 
becomes independent of the mode index n and we can 
reduce the relative correction to the current-current cor- 
relation function and HC to the form: 



was obtained: 



(7 xy 



ni°J(w) 



5T 



(6.13) 



where we have introduced the quantity 

^^^^""^ --iTj2e^^n^[Vom-V,m (6.14) 



gr 



The expression in the RHS of Eq. (|6.14p should be treated 
as a relative correction to the tunneling conductance gT 
of the individual contact due to the tunneling anomaly. 
[The factor 2 in Eq. (|6.13p stands for two contacts ac- 
cording to the square in Eq. (|1.5p for a^J ■] Such 
physical interpretation most clearly arises from the cal- 
culations of interaction corrections to LC. The bare LC 
crlcx (X gr [Eq. (|1.4p ] (in the lowest in gr/go ^ 1 or- 
der) is expressed solely via the tunneling conductance gr 
and the interaction correction to cr^J corresponds to the 
renormalization of gx, since no other physical parame- 
ters are available. The interaction corrections to LC ctS' 
were studied in Refs. [glfiol and at T F the correction 
<5aJ/(c.) to LC 

a.,(c.) = a("i+5aj,^(c.) (6.15) 



gr 



(6.16) 



Therefore, according to the form of Eq. (|6.13p . the cor- 
rection dH^^{uj) indeed corresponds to the renormaliza- 
tion of the tunnelling conductances Gt in Eq. (II. 5p and 
does not affect the Hall resistance Rh of the grain. 



4- Virtual diffusion contribution 511^^ (lo) 

Now let us discuss the correction 611^^ (lu) [Eq. (|6.10p ]. 
The diagram (d), which contributes solely to ^^^'(cl'), 
describes electron diffusion through the central grain. 
The corresponding diffuson D(uj + f7,r,r') enters 
Eq. (|6.10p as the combination 



E 



In 



D ^{uj + fi) - D\_{uj + Q) 



[see Eqs. (j5.4p and (|5.5p ] and therefore contains nonzero 
modes n > Q only. 

As discussed above, one cannot "construct" the prop- 
agator Z)(r,r') [see Eqs. (|5.15l) . (|5.12p ]. which describes 
the real propagation of electron density in metals (basi- 
cally, the classical conduction process) , from the diffuson 
D{u} -\- Vl, r, r') by inserting interaction lines into the cen- 
tral diffuson in the diagram (d): the corresponding dia- 
grams simply cancel each other. This emphasizes virtual 
character of the diffusion: real diffusion is not possible in 
metallic samples, since nonequilibrium electron density 
created in the course of diffusion is screened by Coulomb 
interaction. 

Thus, the correction 51]}^^ {u) describes the process of 
"virtual diffusion" (VD) of electrons through the grain. 



5. Why f2, T ~ Eth are necessary 

We emphasize that it has become possible to identify 
two physically different contributions (5H!^j^(w) [Eq. (|6.9p ] 
and ^n^([j) [Eq. (pTU)) ] to the total correction 6IVxy{uj) 
[Eq. (|6.8p ] only because we have included "high" fre- 
quency range f2 ^ i^xh in our calculations. At frequen- 
cies f2 <C -Exh the expressions in Eqs. (|6.9p and (|6.10p 
acquire the same functional form, since l/(w-f f2-t-7„) w 
1 /7„ = r„ and one cannot distinguish between them. 

Including the frequencies ^ i?Th has complicated the 
calculations significantly. Indeed, for f2 <C i?Th one may 
consider the diagrams (al),(bl),(cl),(d) only. The rest of 
the diagrams are smaller than these ones in fi/Sxh ^ 1 
or uj/Et:\i <C 1 and become comparable to them only at 
^1,0} ^ -ETh- Considering the diagrams (al),(bl),(cl),(d) 
only and neglecting the frequencies w, f2 compared to 
In ^ E-Th in the diffusons l/(w-|-7„) and l/(w-|-f7-t-7„). 
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we would get the total correction SIlxy{i^) [Eq. (|6.8p ] in 
the form 

sn^yH = -2|:^/„T5:0.,of 

n>0 f2>0 

X [3Fo(f^)-2yi(r!)- 1/2(^^)1, (6.17) 

which is just the sum of <5n^jf (w) and m^jf (w) provided 
one neglects u + fl in l/(ti> + f2 + 7„). Clearly, the above 
physical analysis would not be possible based on the form 
of Eq. (ICT71) . 

Another related drawback of considering small fre- 
quencies n <IC i?Th only will be clear in the next sub- 
section, where we turn to the temperature dependence of 
the obtained corrections. 



6. Temperature dependence of the corrections 

Now let us discuss what consequences the corrections 
. (|6.10p have on the Hall conductivity a^y{LLj) 
and, more importantly, on the Hall resistivity 

which is directly measurable experimentally. The con- 
ductivities axy{uj) and a^xi^) are given by Eqs. (|6.1ip 

and (|6.15|) . respectively, and the "bare" HR pl^y by 
Eq. dOiH . 

First of all, since p'^y [Eq. (|5.24p ] simply does not con- 
tain Gt, the tunneling anomaly effect does not influence 
the Hall resistivity. Indeed, according to Eqs. (|6.13p and 
(I6.16p . we obtain that the correction to HR due to TA 
vanishes: 



Pxy 



o xy 



O'xx 



Next, since the analog of VD correction (|6.10p is absent 
for LC axx{^) [see Eq. (I6.15p ] in the leading in gr/go ^ 1 
order^^ , the correction to Hall resistivity due to virtual 
diffusion process is 

Sp^y^ju;) Sa^y^iu;) SU^y^jc.) 



Pxy 



.(0) 



'xy 

Therefore, we obtain that the total correction 

Spxyiu;)^6ply^ic.) + 6pl^^iu;) 
to HR at temperatures T ^ F is due to VD effect only: 

Spxyi^^) _ Scr^y^jc^) 



Pxy 



'xy 



Now, let us discuss the temperature dependence of the 
obtained corrections. The Coulomb potential (|6.ip is 



completely screened and equal to the inverse polariza- 
tion operator. 



2r^ 



for frequencies < gxEc- In the limit of high frequen- 
cies r2 > grEc the Coulomb potential (|6.ip remains un- 
screened and V{n,ci) ~ Ec- The expression 



(6.18) 



in Eqs. (|6.9p and (|6.10p is thus proportional to 1/il for 
57 < grEc- Therefore, the sum of over is logarithmi- 
cally divergent and we have to determine the low and 
high frequency cut-offs. 

In the d.c. limit a; ^ T, the low frequency cut-off for 
the sum over is set by the temperature T (for T 3> F). 
This can be obtained by the analytical continuation of 
Eqs. (|6.9p and (|6.10p to real frequencies cu and taking the 
limit uj <ti T according to 

we do not repeat this standard procedure here (see, e.g., 
Refs. lTlll2^ the integer indices of the Matsubara frequen- 
cies uJn and Qm were written for clarity). 

The high-frequency cut-off is different for TA contribu- 
tion Sli^^ and VD contribution SH"^^ , which is a direct 
consequence of their different physical origin. According 
to Eq. (|6.14p the upper cut-off for 5Il^y is grEc- At 
n > gxEc we have V{^)/n'^ ~ Ec/^'^ and the sum con- 
verges. Therefore, the correction (|6.14p to the tunneling 
conductance gr takes the form^ii^ 



SgriT) 



1 .InP^^" 



gr 



2TTgTd T 



(6.19) 



and for the TA correction to HC from Eqs. (|6.9p and 
(pl^ we get: 



,(0) 
'xy 



1 ,y^^9tEc 



irgxd T 



(6.20) 



The lattice-specific factor 1/d arises as an integral 

1 f a'^d'^q 1 — cos q^a 

C^^y E/3(l - COSqfjo)' 



(6.21) 



For the VD contribution SIlY.^ the summed over 

X y 

f2 expression contains additional f2 in the denomina- 
tor coming from the intragrain nonzero diffusion mode 
1/{lu + fl + 7„). Therefore, the expression is propor- 
tional to 1/n provided not only < gxEc, but also 
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J7 < i?Th- Thus, the upper cut-ofF is the minimum of 
the two quantities, i.e., min ( i^xh j St-S'c)- Calculating the 
sum in Eq. (|6.10p and extracting a^y with the help of 
Eq. (|5T11) . we get 



Sa^y^iT) _ Cd 



In 



min (gT-Ec, -Exh) 



T 



where 



Cd 



a d (\ (1 — cosq2:a)(l — cosg^a) 



5 9/30) 



(6.22) 



(6.23) 



is the lattice form- factor. Appearance of the Thouless 
energy i?Th as an additional cut-off in Eq. (I6.22p reflects 
the diffusive nature of the contribution. Virtual diffusion 
process is suppressed for T > Etu, since in this case the 
intragrain thermal length Lt — y'Do/T < a becomes 
smaller than the size a of the grain. 

We conclude that, the total correction to HR at tem- 
peratures r 3> r is due to virtual diffusion process only 
and equals 



5Pxy{T) 



Cd 



In 



min (ffT-Ec, Et^) 



T 



.(0) J") Airgr 

fJxy ^ xy ^ 

(6.24) 

We see that due to the logarithmic divergence of the 
corrections ()6.20p and (|6.22p one is forced to go to the 
frequencies 51 ^ i?Th in order to get a correct upper cut- 
off, even if one considers the temperatures T <^ i?Th- 
This has direct consequences on physical quantities. The 
upper cut-off also determines the upper bound for the 
temperature range, in which the corrections have the 
InT-dependence of Eqs. ([011)1 and (jO^ . So, HR 
Pxy{T) = p^xy + 5pxy{T) is In T-dcpeudcnt according 
to Eq. ([QI]) for V < T < mmigrEc, Erh), whereas 
LR Pxx{T) = l/a-xx{T) is In T-dependent according to 
Eq. ([grg]) for r < T < gxEc. For T > min(gT-Ec, -Exh) 
the relative correction to HR 6pxy{T)/ p^xy ^ 1/<?t be- 
comes insignificant. The ratio of gxEc and -Exh can be ar- 
bitrary in a real system. In case i?Th QtEc HR Pxy(T) 
is In T-dependent in a narrower range T < T < E^h than 
LR Pxx{T). We emphasize, that perturbative approach 
used by us is applicable as long as the relative correction 
([O^ is small. 

Note that the integrals over q in Eqs. (|6.2ip and 
(|6.23p do not diverge at small qa ^ 1, although Fq = 
2F^^(1 — cosq^a) ^ as — > 0. This is a conse- 
quence of the gauge-invariance mentioned previously. At 
the same time, this is not the case for the density of 
states^. 

The lattice factors 1/d and Cd in 
Eqs. (11201), (IS22), (ESI) are not universal and are 
specific for quadratic/cubic lattice we considered. 
However the In T-dependence itself is a result of the 
screened form of the Coulomb interaction in granular 
metals and is robust to the lattice structure. For a 
different type of lattice the logarithmic dependence of 



Eqs. ((6?20)) . ((O2|) . ((O4)) remains the same, although 
numerical prefactors may be different. 

Moreover, we expect the logarithmic form of 
Eqs. (pT^ . fO^ . lpTM)) to persist even if one allow for 
fluctuations (at least, moderate) of tunneling conduc- 
tances Gt from contact to contact. The reason is that 
the logarithmic contribution arises from the integration 
over frequency: J^^^" dQ/D, = lii{gTEc/T), which is "de- 
coupled" from the integration over quasimomentum q in 
Eqs. (|6.9p . (|6.10p [see Eq. (|6.18p ] in the considered range 
of Q. The tunneling conductance gx enters as an up- 
per cutoff, which precise value with logarithmic accu- 
racy is not important. So, for a realistic array with 
randomly distributed tunneling condcutances the form 
of Eqs. (11201), (IS211), (ESI) should remain, although the 
lattice structure factors 1/d and Cd should be replaced by 
some other distribution-dependent factors of the order of 
unity and gx by some distribution-averaged value. 

We mention in this respect that the logarithmic renor- 
malization of individual conductances for an array with 
randomly distributed tunneling conductances was stud- 
ied by Feigelman, losclcvich, and Skvortsov in Ref. [2^ 
and, indeed, it was shown that the In T-dependence of 
Eq. (|6.19p for the effective tunneling conductance per- 
sists. 



7. Estimate for the contribution from nonzero interaction 
modes. 



As we claimed above, coordinate-dependent intra- 
grain interaction modes Sijv{il,ri,r'^) [Eq. (|4.16p ] of 
the screened potential V{^l,ri,r'j) [Eq. (|4.19p ] give a 
smaller contribution than the zero-mode part V{fl,i,]) 
[Eq. (11201)] in the relevant range fl,T < ETh- We pro- 
vide an estimate here. 

Let us revise, for example, the diagrams (al),(a2) in 
Figs.fT^andfinitaking now nonzero modes of the screened 
interaction into account. For the block in the right grain, 
we have 

J drdr'D{n,s,r)V{n,r,r')D{n,r',s) = 



n>o 



O>0 



n>0 



(fl-|-7„)2 



The first term in the RHS is the contribution from the 
zero-mode interaction Vo{^l) we had before, whereas the 
second one is the contribution from nonzero modes u„(r2), 
which we want to estimate now (s is a point on the con- 
tact). Since Vb(ri) ~ Sfl/T for Q < grEc, the first term 
gives the logarithmic contribution (j6.19p to the conduc- 
tivity: 



1 

gr 



In 



T 



(6.25) 
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As Vn{n) ~ (fi + 7„)/7„ for n < niin(l/ro, afj = I/tq, 
the contribution from nonzero modes is estimated as 



1 



1 



TTsyT, 



1 



I'D, 



dn 



Jt 



dq 



iPFiy 



[1 - Q;v/Tomax(r, E'xh)] 



(6.26) 



(here a ^ 1 is a cutoff-dependent number). The ob- 
tained contribution from nonzero modes is nothing else 
but the tunnehng anomaly in 3D caseiii^^i^ (since the 
grains are three-dimensional) , the square root describing 
well-known singularity. We see, that the contribution 
(|6.26|) from nonzero modes is smaller than the contribu- 
tion (I6.25|) from zero modes if the condition 



90 /gr > a/ 1 



(6.27) 



is met. Since go » gr [Eq. (|1.3p ]. for ballistic grains 
{I ~ a) this is always the case. More important, however, 
is that in the relevant temperature range T < i^xh the 
contribution (|6.26p is T-independent. Therefore, even 
if the condition (|6.27p is not well met, nonzero modes 
give an inessential T-independent renormalization of the 
bare quantities for T < i?Th not affecting the overall T- 
dependence of the HC and HR in that range. 



B. 



"Low" temperatures T < F 



Now we want to include the region of temperatures 
T < r of the order or smaller than the escape rate T into 
our considerations. In this regime the thermal length for 
the intergrain motion = y'Ta^/T ^ a is greater than 
the size of the grain a and quantum phenomena can occur 
not only inside the grains, but also at spatial scales much 
exceeding the grain size. 

The technical complication in considering T < F in- 
stead of T ^ F is that we have to take tunneling fully 
into account. At the same time we still have a small 
parameter F/ii'Th ^ 1 [Eq- (|1.2p ]. which allows us to 
neglect nonzero intragrain modes compared to the zero 
modes as long as we consider the frequencies fl <^ -ETh- 

As we have discussed in the previous section, even for 
temperatures T <C ETh considering ^ i?Tii in the ex- 
pression 



n>o 

for the correction to HC is necessary in order to obtain 
a correct upper cut-off for logarithmically the diverging 
quantities. At the same time for ~ E'xh the results 
must match with those of the previous section, since for 
^ F one can neglect tunneling in the expressions. 
Therefore, to simplify our calculations we do assume 
<C ETh in this section. Having obtained the results 



limited by <C -E^Th, we determine the upper cut-off by 
matching them with the results of the previous section 
in the range T ^ ft ^ i?Th, where both results are ap- 
plicable. 

The diagrams for T <r can be obtained from the di- 
agrams in Figs. [T^ and [T51 for T ^ F by including higher 
orders in tunneling. As we restrict ourselves to low fre- 
quencies il <C ETh, only four diagrams (al),(bl),(cl),(d) 
should be considered and the rest of the diagrams are 
smaller in f2/i?Th ^ 1- Taking tunneling into account 
we make sure that in each diagram only one grain con- 
tains nonzero dijfusion modes necessary to have a nonva- 
nishing contribution to the Hall current, whereas in the 
rest of the grains only the zero mode is retained. The 
zero-mode interaction potential V{Q, i, j) is now taken in 
the form (|4.20p . Considering the range llj,^! ETh, we 
also neglect the frequency dependence of nonzero modes 
now: l/(w + + 7„) = l/7„ and 1/ (w + 7„) = l/7„. 



1. Short-scale contribution 

Including higher orders in tunneling in diagrams (al), 
(bl), and (cl) in Fig. [T^ is straightforward. One has 
to substitute the zero- mode diffusons Sij/^ renormal- 
izing the interaction vertices by their form I?o(ri,i,j) 
[Eq. (|4.18p ]. which takes tunneling into account: 



1 



fi + Fo 



This results in the replacement of V{fl,i,j)/ft^ in the 
expressions for Xif'iuj^n), i = al,bl,cl, (see Table HH by 

k,l 



Vin,q)=VliQ,q)V{n,q). 



Consequently, instead of Eqs. (|6.7ap . (|6.7bp . and (|6.7cp 
we get [Fig.[lla),(b),(c)] 



aaA^(w,r!) = -4T^o(f^)- 



In 



(6.28a) 
(6.28b) 



a,Xl{iu,n)^~2Vo{n) — , (6.28c) 

In 

where Voifl) = V{n,i,i) and Vi{fl) = V{fl,i + e^,i) ^ 

The diagram (d) in Fig. [T^ has to be considered care- 
fully. It contains three diffusons: one "central" diffu- 
son l/(w -|- O -|- 7„) describing diffusion through the sin- 
gle grain and two "adjacent" diffusons 1/ft renormaliz- 
ing the interaction vertices. In the general case of ar- 
bitrary compared to F temperatures T, the diagram (d) 
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FIG. 14: Diagrams for the interaction corrections to Hall con- 
ductivity arising from spatial scales of the order of the grain 
size ("short-scale" contributions) at arbitrary compared to 
r temperatures T. Gray blocks denote nonzero-mode intra- 
grain diffusons D{Q,r,v') [Eq. (I4.9|l ]. whereas rendered with 
lines blocks denote zero-mode diffusons 'DoiO., i, j) [Eq. (|4.18|l ] 
of the whole system. 



for the current-current correlation function nab(w,i,j) 
corresponds to the process of virtual diffusion, when an 
electron, "created" at the contact (j + b, j) by the applied 
bias, gets diffusively, without additional applied bias, to 
the contact (i -I- a, i), thus contributing to the current. 
Only in the limit T 3> F the main contribution comes 
from the closest contacts of a single grain. This is a vir- 
tual process, since for a real electron its charge would be 
screened. 

Accounting for tunneling in the diagram (d) in Fig. 1121 
one should, in principle, replace each one of the three dif- 
fusons by the "exact" diffuson V [Eq. (|4.2p ]. However, 
in the lowest nonvanishing in gT/9o oc T/E^h <^ 1 or- 
der it is sufficient to leave the non-zero-mode diffuson 
D{0, r, r') [Eq. (|4.9[) ] in only one of the grains and retain 
only the zero modes l/(|i7|V) in the rest of the grains. 

One possibility is to leave the "central" diffuson as it is 
in Fig. [121 i-e., as a non-zero-mode diffuson D{llj, r, r') of 
a single grain, and to take "adjacent" diffusons as zero- 
mode diffusons 'Do{fl,i,j) [Fig. [Hf d)]. This gives the 
contribution analogous to Eq. (|6.7dp : 



adXt{Lo,n) = 2^2 (r!) 



1 

In'' 



(6.28d) 



where V2(ri) = V{il., i + e^;, i + ey). The rest of the di- 
agrams arising from the diagram (d) in Fig. [T^] are dis- 
cussed in the next section. 

Le t us su m the contributions (|6.28a|) . (j6.28bp . (|6.28cp . 
and (|6.28dp . For the corresponding correction 



7l>0 



O>0 



a,b,c,d 



to the bare current-current correlation function H^xyi^^) 



[Eq. ((5TT|) ] we obtain 



n>0 



1 



X —[Wo{n)^2Vi{n)-V2{n)].{6.29) 

In 



Compare the correction (|6.29p with the correction (|6.17p 
obtained in the regime T ^ F, in which fJ, T <C -Exh was 
also assumed. Equation (|6.29[) differs from Eq. (|6.17[) 
only by having 'Do{n,q) = + Fq) instead of l/ft 

for the diffusons rcnormalizing the interaction vertices. 
This difference is relevant at T, f2 < F, when tunneling 
comes into play. Nevertheless, the correction (|6.29p still 
arises from the spatial scales of the order of the grain size 
a, even if T <C F. At T, ^ F, when one can neglect 
tunnehng in Eq. (|6.29p [I?o(f^, q) 1/^]: the corrections 
and (|07| coincide. 
To get a temperature dependence of the correction 
(I6.29P we note that the summed over Q expression is 
still proportional to in the wide range F < O < 
grEc, giving a logarithmic contribution. The question 
about the lower cut-off is easily resolved: it equals T 
at r > F and F at T < F. Thus, the lower cut-off 
is max(T, F). To get a correct upper cut-off we separate 
Eq. (j6.29p artificially into two parts according to the form 

of Eqs. mm. mjm : 



where 



and 



,5t 



n>0 



O>0 



1 1 



(6.30) 



In 



n>0 n>o 
X ^[v2{n) + Vo{n)-2ViiQ) .(6.31) 



Keeping in mind that Eqs. (|6.30p and (|6.3ip must match 
with Eqs. (|6.9p and (|6.10p a.t fl ^ -E^Th, we have to at- 
tribute a cut-off grEc to SH'^^{uj) and mm{gTEc, ETh) 
to SII^^{uj). Doing so, for the corresponding corrections 
to HC we obtain 



and 



,(0) 



irgrd 



In 



grEc 



max(T, F) 



for T < grEc 

(6.32) 



Cd 



47r5T 

for T < min((/T iJc, i?Th)- 



In 



min(5T£'c,£^Th) 



max(r, F) 



(6.33) 
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Let us write the total correction 5axy to HC 



V, e.g., that Voin^i - j) = Voifl^ix - jx 



jy) 



7 - (t(°) 
'xy — <-'xy 



'xy 



as 



Saxy = 5(7, 



TA 



5a 



VD 



5C7. 



AA 



(6.34) 



where iJcrJ^^ and 5ax^ are given by Eqs. I|6.32p and (|6.33p 
and 5a x^ comes from the rest of the diagrams, arising 
from the diagram (d) in Fig. [T^ when we take tunneUng 
into account. We consider the latter "large-scale" contri- 
bution now. 



2. Large-scale contribution 

As explained above, in order to generalize the diagram 
(d) in Fig. [12] to include the temperatures T < T one 
has to replace all three difFusons in it by the zero-mode 
diffusons I?Oi a-nd then insert a non-zero- mode part D 
[Eq. gH)] of the diffuson D [Eq. gT])] into one of the 
grains. Having considered the diagram (d) in Fig. 1141 we 
are left with two following possibilities. 

(i) One can insert the non-zero-mode part D [Eq. (|4.9p ] 
somewhere into the middle of one of the three difFusons 
Vq. One can show straightforwardly that all the contri- 
butions from such diagrams, when summed up, cancel 
each other exactly. 

(ii) The less trivial possibility is to consider nonzero 
modes D in the grains directly adjacent to the contacts 
{1 + ex, i) and (j -|- e^, j), corresponding to external cur- 
rent vertices for the correlation function \ixy{i^, i, j), i.e., 
in the grains i -|- e^, or i for contact (i -I- e^, , i) and in the 
grains j + or j for contact (j -I- e^, j). However, the 
total contribution from such diagrams also vanishes: for 
each diagram of this type there exists another diagram, 
the contribution of which is exactly opposite and, con- 
sequently, their sum is zero. One of such pairs are the 
diagrams (dl) and (d2) shown Fig. [151 their contributions 
to the current-current correlation function being 

9I. 



X J2 r^o(t-' + ^, i, j + ex)V{n,i - ey,i + e^), 



5Kico) 



9t 



ri?o(w + fl, i, j - ey)V{n,i - e:r, i + e^;). 



Here D> is the non-zero-mode intragrain diffuson con- 
necting contacts in the counterclockwise direction, i.e., 
in the direction of the edge drift. Using the translational 
invariance and symmetry in each component of the zero- 
mode diffuson Vq and the screened Coulomb interaction 



, jx — ixjiy — jy), we obtain that the sum over j in 
5IIxy{i^) is identical to the one in 5Ilxl{u)), and thus the 
contributions cancel each other: 



5U%iu;) + 5Iiry{co) = 0. 

Therefore, we obtain that the total contribution from 
the diagrams of types (i) and (ii) vanishes identically: 



5a. 



AA 



0. 



(6.35) 



As a result, all nonvanishing contributions to HC at 
arbitrary compared to F temperatures T are given by 
the diagrams (a),(b),(c),(d) in Fig. [TH which lead to the 
corrections (|6.32p and (|6.33p . Note that these contribu- 
tions arise from the spatial scales of the order of the grain 
size a, even for temperatures T <C F. On the contrary, 
the eventually vanishing contributions of the diagrams 
of types (i) and (ii) [as (dl) and (d2) in Fig. [T3] arise 
from the spatial scales exceeding a. The reason is that in 
these diagrams the contacts with external tunneling ver- 
tices are connected by the diffuson 'Dq{Q,, q), which "size" 
in the re al spac e is determined by the thermal length 
= yTa?/r. In case of low temperatures T ^ F, 
the thermal length ^ a can exceed the grain size a 
significantly. 

It is always instructive to compare the results for a 
granular metal with those for an ordinary homogeneously 
disordered metal (HDM). For quantities arising from 
large spatial scales (i.e., much greater than the grain size 
a for granular metals and the mean free path for dif- 
fusive HDMs) one expects correspondence between the 
two, since at such scales the microscopic structure of the 
system becomes irrelevant. The first-order Coulomb in- 
teraction correction 5axy to Hall conductivity of HDM 
was first studied in Ref. ^ and the correction was found 
to vanish, 5axy = 0. So, indeed, our result (I6.35|) for 
the "large-scale" contribution 5a^^ agrees with that for 
HDMs. 

Note that even our approach of calculating 5a^,^ is 
quite similar to that of Ref. d. The authors of Ref. [^cal- 
culated 5axy perturbatively in magnetic field H (assum- 
ing ujhtq ^ 1) by inserting "magnetic vertex" — ^Ap 
in all possible ways into the diagrams for zero magnetic 
field, (i) Insertions of magnetic vertex into the difFu- 
sons were found to cancel and (ii) insertions of magnetic 
vertex into the block of Green functions at the current 
vertices (Fig. 6 in Ref. 0) were found to cancel. These 
two steps resemble those (i) and (ii) of our approach, 
insertion of magnetic vertex corresponding to insertion 
of non-zero-mode intragrain diffuson D, which contains 
all information about magnetic field, into the zero-mode 
"large-scale" diffusons Vq. 

We want to stress that the "large-scale" AA contribu- 
tion 5a^^ is not any different physically from the "short- 
scale" VD contribution 5a^^ . They both correspond to 
the process, when electron gets diffusively from the con- 
tact (j -\- By,'}) to the contact (i -I- e^,, i). It just happens. 
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FIG. 15: Diagrams for the interaction corrections to Hall conductivity arising from large spatial scales at arbitrary compared 
to r temperatures T. The diagrams (dl) and (d2) cancel each other leading to the vanishing "Altshuler-Aronov" correction to 
the Hall conductivity: 5a^^ — 0. 



that the contribution Sa^^ from the diffusion processes 
through a single grain does contribute to HC, whereas 
the total contribution 6a^^ to HC from the diffusion 
through more than one grain vanishes. 

As a result of Sees. I VI B II and I VI B 2l we obtain that 
the total correction Sa^y [Eq. (|6.34p ] to HC at arbitrary 
temperatures T is given by two short-scale contributions 
foj/ [Eq. dOa ] and Sa^y^ [Eq. (I05D ]. whereas the 
large-scale contribution Sa^^ = [(Eq. (|6.35p ] vanishes 
in agreement with the theory of HDMs. For the discus- 
sion of the corresponding corrections 

to the Hall resistivity 

"'^y I'D) I s- 

P^v = -;jr = Pxv +^P^y 

XX 

we refer the reader to the Results section |ll] starting from 
Eq. (lO) . 

VII. CONCLUSION 

In conclusion, we presented a theory of the Hall ef- 
fect in granular metals. In spite of its importance this 
question has not been addressed before. It turned out 
that considering only zero intragrain spatial harmonics 
that was very successful in describing the longitudinal 
conductivity^ii^ is not sufficient for calculating Hall con- 
ductivity and we were forced to take nonzero harmonics 
into account. Proceeding in this way we have shown that 
at high enough temperatures the Hall resistivity is given 
by the classical expression, from which one can extract 
the effective carrier density of the system. At lower tem- 
peratures, quantum effects come into play, their most sig- 
nificant effect being logarithmic temperature-dependent 
contributions to the Hall conductivity and resistivity due 
to Coulomb interaction, which are absent in ordinary dis- 
ordered metals. 



We emphasize, however, that Eqs. (|2.4p . (|2.5p give the 
first-order interaction corrections and the result is valid 
as long as these correction are small. One could try to 
account for higher orders using a renormalization group 
analysis (since the contributions are logarithmic), for ex- 
ample, like the one used for the longitudinal conductivity 
in Ref. iQj. In order to write down proper renormalization 
group equations, methods more sophisticated than the 
present diagrammatic approach are needed. 

Concerning the experimental situation related to 
the theory developed here, the logarithmic dependence 
Pxx{T) = i?i — i?2lnT of the longitudinal resistivity 
in good conducting {qt ^ 1) granular materials, cor- 
responding to the logarithmic renormalization^ii^ of the 
integrain tunneling conductance gr [ Eg. (|6.19p ]. has been 
observed experimentally (see Refs. |27| ). Clearly, mea- 
surements of the Hall resistivity pxy of such granular 
samples could be also performed and our theory could 
be thus tested. Unfortunately, the known to us exper- 
imental papers (Refs. [H) on convetional Hall effect in 
granular materials mostly deal with the systems in the 
regime of low tunneling conductance <C 1, opposite 
to the metallic regime ^ 1 studied by us, which does 
not allow us to make a detailed comparison now. We 
mention that our theory may be also applied to indium 
tin oxide (ITO) materials (see, e.g., Refs<^). Another re- 
lated effect is the anomalous Hall effect in ferromagnetic 
granular materials2£. 

We hope that more experiments on this subject will be 
done in the nearest future and that Hall measurements 
will evolve into an important method of characterization 
of granular materials. 
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FIG. 16: Diagrams for the current- (radius vector) correlator 
{jarp) [Eq. (|A6|l ]. if one directly expands Green functions to 
the linear order in vector potential A. Fermionic lines denote 
the Green function [G{e,'p)]^^ = ie — ^(p) + j^sgne of a 
bulk metal with H = 0. (a) Magnetic-field-independent part 
of (jarp) giving the LHS of the boundary condition (|4.6p . (b) 
Linear in magnetic field part of {jar/}) obtained by insert- 
ing "magnetic vertex" — ^Ap in all possible ways into the 
diagram (a) and giving the RHS of Eq. (|4.6p . 



APPENDIX A: BOUNDARY CONDITION 

To obtain the boundary condition for the difFuson 
D{llj, r,r') [Eq. (|4.ip ]. we recall its physical meaning: In 
the real-time representation, the quantity 

/ + 00 
c^^I'e-"^*[DK,r,r')]|^c.„^.D+^o, 
-oo w„>0 

gives the probability density to find an electron at point 
r at time t provided it was at point r' at time t = 0. 
Therefore, according to the formal definition (|4.ip of the 
diffuson, we can write down the probability current cor- 
responding to the diffusion process as 

r, r') = -^(jr [Gie + r, r')g{e, r', r)])^;, {e+Lo)e < 0, 
Zni' 

(Al) 

where 



jr[e(£ + c.,r,r')^(e,r',r)] = —[^(e, r', r)(-zVr)^(£ + c., r, r') + 

+ gie + Lu,r,r'){iS/,)g{e,r',r)]-—A{r)g{s+LO,r,r')g{e,r',r) (A2) 

mc 

I 



is the current operator acting on the product of two 
Green functions and A(r) is a vector potential corre- 
sponding to magnetic field H. Since electron cannot es- 
cape from an isolated grain, the normal component of the 
current j = r, r') must vanish at the grain boundary: 

(n,j)|rGS = 0. (A3) 

Here, n is a unit vector normal to the grain boundary and 
pointing outside the grain. Equation (jASj together with 
Eqs. and ([X2|l gives a general form of the boundary 
condition for the diffusion propagator D(w,r,r'). Fur- 
ther simplifications depend on the model used. In the 
case of white-noise disorder, using the integral equation 
(14. 4p for the diffuson, we obtain 

j(w,r,r') = jr[i?o(^,r,r')] 

+ - / dxjr[i?o(^,r,x)]i?(cc;,x,r')(A4) 

To J 

where jr acts on the Green functions of Dq{llj, r,r') = 
G(e-hw,r,r')G(e,r',r) as in Eq. (IMl) . 

We now exploit the diffusive hmit. Since Do{Lj,r,r') 
varies on the spatial scale Z <C a, we can (i) neglect the 
first term in the RHS of Eq. (|A4|) : (ii) expand £>(tj, x, r') 
writing it as 

D{lj, X, r') » D{lu, r, r') + (x - r)p^,pD{u;, r, r') 



to get 

ja(w,r,r') ^ — (jar^)Vr/3l?(w,r,r'), (A5) 
To 

where 

{iarp) ^ J dxjra[Do{Lj,r,x)]{x- r)f3 (A6) 

is the current- (radius vector) correlation function, (iii) In 
the diffusive limit, details of the boundary scattering be- 
come unimportant as soon as we move away from the 
boundary into the bulk of the grain over the distance of 
the order of the mean free path I. Since D{oj,r,r') in 
Eq. (jA5P varies on the scales a ^ I, the condition (jASP 
may be evaluated not exactly at the boundary, but at 
some point a few I away from it. This allows us to use 
for Green functions G in Eq. (|A6p their expressions for 
the bulk. 

The correlation function {jaTfj) can be calculated via 
diagrammatic technique either by directly expanding 
Green functions in vector potential A(r) (Fig. [TB)) or us- 
ing an explicitly gauge-invariant approach developed by 
Khodas and Finkel'stein in Ref. [ij. Proceeding in either 
way, in the linear order in H we get 

Oaf/3) = C {Sap + "^0/37-^7) ' (^'^) 
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where eap-y is the totally antisymmetric tensor, txyz = 1, which can be easily expressed in the form of Eq. (|4.6p . 
and C = —{21: /2>)vP is an irrelevant for the boundary 
condition (|A3P prefactor. Inserting Eqs. (|A5P and (jA7P 
into Eq. (jASp . we obtain the boundary condition for the 
diffuson D = D{uj, r, r') in the presence of magnetic field: 




^ B.L. Altshuler and A.G. Aronov, in Electron- Electron In- 
teractions in Disordered Systems, edited by A.L. Efros and 
M. Pollak (Elsevier, Amsterdam, 1985). 

^ P.A. Lee and T.V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

^ I.S. Beloborodov, K.B. Efetov, A.V. Lopatin, and V.M. 
Vinokur, Rev. Mod. Phys. 79, 469 (2007). 

* M.Yu. Kharitonov, K.B. Efetov, Phys. Rev. Lett, accepted, 
,cond-mat/0609736 

^ Yu.V. Nazarov, Sov. Phys. JETP 68, 561 (1989). 

® L.S. Levitov, A. V. Shytov, JETP Lett. 66, 214 (1997). 

I.S. Beloborodov, K.B. Efetov, A.Altland, and F.W.J. 

Hekking, Phys. Rev. B 63, 115109 (2001). 

* B.L. Altshuler, D.E. Khmelnitski, A.I. Larkin, and P.A. 
Lee, Phys. Rev. B 22, 5142 (1980). 

^ K.B. Efetov and A. Tschersich, Europhys. Lett. 59, 114, 

(2002) ; Phys. Rev. B 67, 174205 (2003). 

^° I. S. Beloborodov, K. B. Efetov, A. V. Lopatin, and V. M. 
Vinokur, Phys. Rev. Lett. 91, 246801 (2003). 

B. L. Altshuler and A.G. Aronov, Sov. Phys. JETP 50, 968 
(1979). 

I.S. Beloborodov, A.V. Lopatin, and V.M. Vinokur, Phys. 
Rev. B 70, 205120 (2004). 

Ya.M. Blanter, V.M. Vinokur, and L.I. Glazman Phys. 
Rev. B 73, 165322 (2006). 

C. Biagini, T. Caneva, V. Tognetti, and A. A. Varlamov, 
Phys. Rev. B 72, 041102(R) (2005). 

M.Yu. Kharitonov, K.B. Efetov, in preparation. 
1^ H. Fukuyama, J. Phys. Soc. Japan 49, 644 (1980). 
" M. Khodas, A. M. Finkel'stein, Phys. Rev. B 68, 155114 

(2003) . 

A. A. Abrikosov, L.P. Gorkov, and I.E. Dzyaloshinski, 
Methods of Quantum Field Theory in Statistical Physics 
(Prentice-Hall, Englewood Cliffs, New Jersey, 1963). 

I.L. Aleiner, P.W. Brouwer, and L.I. Glazman, Phys. Rep. 
358, 309 (2002). 

L.D. Landau and E.M. Lifshitz, Electrodynamics of Con- 
tinuous Media, Course of Theoretical Physics, Vol. 8 (Perg- 
amon. New York, 1975). 
^1 A.M. Finkelstein, Sov. Phys. JETP 57, 97 (1983); C. 
Castellani, C. Di Castro, P.A. Lee, and M. Ma, Phys. Rev. 
B 30, 527 (1984). 

G. Zala, B.N. Narozhny, and I.L. Aleiner, Phys. Rev. B 
64, 214204 (2001). 

L.P. Gorkov, A.I. Larkin, and D.E. Khmelnitskii, Sov. 
Phys. JETP Lett. 30, 228 (1979); S. Hikami, Phys. Rev. 
B 24, 2671 (1981). 

M.V. Feigel'man, A.S. loselevich, M.A. Skvortsov, Phys. 
Rev. Lett. 93, 136403 (2004). 

B. L. Altshuler and A.G. Aronov, Solid State Commun. 30, 
115 (1979). 

B.L. Altshuler and A.G. Aronov, Solid State Commun. 46, 



429 (1983). 

R.W. Simon, B.J. Dalrymple, D.V. Vechten, W.W. Fuller, 
and S.A. Wolf Phys. Rev. B 36, 1962 (1987); 
H. Fujimori, S. Mitani, S. Ohnuma, T. Ikeda, T. Shima, 
and T. Masumoto, Mater. Sci. Eng., A 181, 897 (1994); 
A. Gerber, A. Milner, G. Deutscher, M. Karpovsky, and 

A. Gladkikh, Phys. Rev. Lett. 78, 4277 (1997); 

L. Rotkina, S. Oh, J. N. Eckstein, and S. V. Rotkin Phys. 
Rev. B 72, 233407 (2005). 
^® X.X. Zhang, Chuncheng Wan, H. Liu, Z.Q. Li, Ping Sheng, 
J.J. Lin, Phys. Rev. Lett. 86, 5562 (2001); 

B. Bandyopadhyay, P. Lindenfeld, W.L. McLean, and H.K. 
Sin Phys. Rev. B 26, 3476 (1982); 

M. Nissim and R.L. Rosenbaum Phys. Rev. B 40, 10629 
(1989); 

T.T.M. Palstra, R.C. Haddon, A.F. Hebard, and J. Zaanen 
Phys. Rev. Lett. 68, 1054 (1992); 

X.N. Jing, N. Wang, A.B. Pakhomov, K.K. Fung, and X. 

Yan Phys. Rev. B 53, 14032 (1996). 
'^^ J. Ederth et al, Phys. Rev. B 68, 155410 (2003); J. Ederth 

et al, Thin Solid Films. 445, 199 (2003). 
^° P. Mitra, A. F. Hebard, K.A. Muttahb, P. Woelfle, 

cond-mat/0606215 

Note that although the granular array may be two- {d = 2) 
or three-dimensional (d = 3), the grains themselves are 
three-dimensional, and n is a three-dimensional density. 
Dealing with the Hall transport, we do not discuss one- 
dimensional case of granular "wires" in this paper, for 
which the "Altshuler- Aronov" correction 5a^^ is also di- 
vergent. 

Roughly, the reason is that diagram for the bare LC a^J 
[Eq. (|1.4|l . Fig. [3jd)] does not contain diffusons, contrary 
to the diagrams for HC aSj (Figs. I6I8P , and consequently, 
the diagrams for the Coulomb interaction corrections de- 
scribing the "virtual diffusion" process do not arise. See 
also the footnote on page 1211 

Although we term the propagator D(tj,r,r') defined by 
Eq. (|4.H) as "diffuson" , there is yet no need to assume the 
diffusive limit {I <^ a) for these general considerations and 
they are valid in the ballistic case {I > a) as well. 
To avoid misunderstanding, we emphasize that zero modes 
drop out from the expression for the classical conductivity 
only for the system with identical tunneling conductances 
Gt of all contacts. If conductances are not equal, then, 
for example for LC in the limit Gt/Gq ~ 0, one needs 
to consider not only the contribution from a single given 
contact [a simple diagram (d) in Fig. but also those 
from all other contacts, which have to be connected to a 
given contact by zero-mode diffusons. Briefly, zero modes 
take care of the fact that conductances of different contacts 
are not equal, while nonzero modes take care of the finite 
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resistances of the grains themselves (that Gt/Gq 7^ or 
GtRh 7^ 0). Our aim is to discuss the latter point and to 
show that this is crucial for the Hall effect. 
In fact, the Coulomb interaction has to be also taken into 
account in a certain way to get a correct classical expres- 
sion for (Tab(<i') at nonzero uj. This point will be discussed 
in detail in the case of Hall conductivity in SecFVl 
We neglect the dispersion a^^iuj) at ~ 1/ro of the grain 



conductivity itself in these considerations. 
In fact, the correction 5a^^ analogous to Sa]^^ , also exists 
for LC. However, it is a correction to the term —Gt/Gq 
in Eq. (|4.3p but not to Gt- Therefore it contains an addi- 
tional smallness Qt/qo ^ 1 compared to Jcrjjf and can be 
neglected: Sa^^/Sa^^ oc gr/go < 1- 



